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SECTION  1 

INTRODUCTION 


A  recent  paper  [1]  considered  various  trade-offs  in  the  use  of  spread  spectrum 
communications  over  the  wide  band  (e.g.,  1-MHz)  HF  channel.  It  was  concluded  that 
direct  sequence  pseudo-noise  (DSPN)  spreading,  with  interference  excision  and  adaptive 
matched  filtering  at  the  receiver,  can  provide  an  efficient  communication  technique,  with 
links  in  the  1000-3000  km  range  being  particularly  favorable. 

This  report  examines  the  effectiveness  of  four  different  Rake  techniques  as  a  means 
to  achieve  adaptive  matched  filtering  of  DSPN  signals  received  over  the  wide  band  HF 
channel.  The  Rake  technique  was  conceived  by  R.  Price  and  P.  Green  over  thirty  years 
ago  [2]  as  a  means  of  combating  multipath  and  was,  in  fact,  built  and  tested  over  a 
narrowband  (10  kHz)  HF  channel.  Current  technology  makes  viable  the  application  of 
the  Rake  technique  to  the  wide  band  HF  channel. 

The  presence  of  strong  narrowband  interfering  signals  throughout  the  HF  band 
requires  that  any  DSPN  wide  band  HF  communication  system  employ  a  narrowband 
interference  excisor  preceding  the  data  demodulation.  Based  upon  actual  interference 
measurements,  it  is  shown  in  [3]  that  only  a  relatively  small  percentage  of  a  1-MHz  band 
may  need  to  be  excised  to  achieve  large  improvements  in  signal-to-noise  ratio  (SNR). 
As  long  as  the  percentage  of  the  band  excised  is  small,  the  effect  of  the  excisor  on  the 
DSPN  signal  is  a  small  reduction  in  amplitude  plus  the  creation  of  a  small  amount  of 
self  noise.  However,  this  self  noise  will  be  negligible  compared  to  the  noise  at  the  Rake 
combiner  output.  The  performance  analysis  presented  here  does  not  include  the  effect 
of  the  excisor  per  se,  since  it  expresses  performance  in  terms  of  the  value  of  Eb/N0 
(energy/bit/noise  power  density)  at  the  input  to  the  Rake  processor.  Any  self  noise 
due  to  the  excisor  is  neglected  and,  in  addition,  the  residual  noise  at  the  excisor  output 
is  assumed  white.  The  validity  of  this  latter  assumption  has  not  yet  been  established, 
although  measurements  at  The  MITRE  Corporation,  Bedford,  MA  facility  do  show  an 
essentially  white  noise  behavior  at  the  excisor  output. 

Prediction  of  modem  performance  over  the  HF  channel,  wide  band  or  not,  is  a 
difficult  task  due  to  the  random  and  unpredictable  nature  of  various  HF  channel  pa¬ 
rameters  and  to  the  lack  of  definitive  statistical  characterization  of  these  parameters. 
As  discussed  further  in  section  2,  we  sidestep  these  issues  by  confining  our  performance 
calculations  to  non-overlapping,  purely  dispersive  propagation  modes  and  conditioning 
our  calculations  on  specified  mode  strengths  and  Doppler  shifts.  Nonetheless,  these  re¬ 
sults  are  useful  in  providing  comparative  modem  performance  results.  In  addition,  once 
an  adequate  statistical  description  of  mode  strengths  and  Doppler  shifts  is  available, 
the  conditional  results  presented  here  may  be  unconditioned  by  averaging  over  these 
parameters. 
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Section  2  summarizes  the  wide  band  HF  model  used  in  the  performance  analysis 
and  section  3  briefly  describes  the  four  Rake  modems  analyzed.  Some  performance  re¬ 
sults  are  presented  in  section  4.  The  detailed  analysis  supporting  the  results  is  presented 
in  the  appendix. 
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SECTION  2 

CHANNEL  MODEL 


The  HF  propagation  medium  can  be  represented  as  a  linear,  randomly  time- variant 
operation.  From  innumerable  measurements  and  propagation  analyses,  it  is  known  that 
the  received  HF  signal  can  be  expressed  closely  as  the  sum  of  contributions  from  one  or 
more  propagation  modes  (e.g.,  see  [4]).  Each  propagation  mode  itself  may  be  regarded 
as  a  time-variant  linear  operation.  This  linear  operation  will  have  different  statistical 
structures,  depending  upon  the  ionospheric  conditions  in  the  corresponding  layer.  Thus, 
if  the  percentage  change  in  electron  density  is  small  enough,  the  propagation  is  called 
quiet  or  normal,  and  a  mode  acts  over  a  MHz  bandwidth  approximately  like  an  all-pass 
filter  with  a  somewhat  nonlinear  group  delay  vs.  frequency,  a  slowly-changing  mean  time 
delay,  complex  amplitude,  and  Doppler  shift.  For  path  lengths  in  the  1000-3000  km 
range,  and  operations  somewhat  below  the  maximum  useable  frequency  (MUF),  the 
group  delay  is  nearly  linear  vs.  frequency  and  dispersions  are  usually  not  greater  than 
100  /xs/MHz.  Although  Doppler  shifts  can  be  of  the  order  of  1  Hz,  this  usually  occurs 
at  dawn  or  dusk.  For  shorter  path  lengths  these  numbers  increase.  When  the  electron 
density  fluctuations  become  sufficiently  large,  the  mode  acts  like  a  scatter-type  channel 
with  an  associated  Doppler  spread  and  a  multipath  spread  caused  by  random  fluctu¬ 
ations  of  received  energy  over  a  range  of  path  delays.  The  mode  may  then  be  called 
disturbed.  Doppler  spreads  of  tens  of  Hz  and  multipath  spreads  of  tens  of  milliseconds 
may  occur  in  extreme  cases.  Such  variations  in  mode  properties  with  the  degree  of 
electron  density  fluctuations  have  been  discussed  by  Booker  et  al  [5] [6]  for  the  F  layer. 

For  the  purpose  of  our  analysis,  the  channel  per  se  must  be  regarded  as  a  com¬ 
posite  filter  consisting  of  the  cascade  of  all  transmitter  filtering  operations,  the  HF 
propagation  medium,  and  all  receiver  filtering  operations.  This  view  makes  clear  that 
the  ability  to  resolve  propagation  modes  depends  not  only  on  how  close  these  modes 
are  in  propagation  delay,  but  also  on  the  impulse  response  duration  of  the  transmitter 
and  receiver  filters.  In  evaluating  the  performance  of  the  Rake  modem,  it  is  useful  to 
place  the  wide  band  HF  channel  into  one  of  the  following  categories: 

(a)  non-overlapping,  purely  dispersive  propagation  modes 

(b)  partially  overlapping,  purely  dispersive  propagation  modes 

(c)  scatter  type  or  combination  scat  ter /dispersive  propagation  modes. 

Case  (a)  is  the  category  that  corresponds,  ideally,  to  the  use  of  wide  band  HF 
under  normal,  i.e.,  non-disturbed  HF  propagation  conditions.  All  modes  are  resolvable 
and  have  complex  amplitudes  that  vary  very  slowly  and  do  not  fade  deeply  (e.g.,  time 
constants  exceeding  10  seconds  and  CL  fCW  il  B  fade  depth).  However,  Doppler  shifts  exist 
with  values  that  may  be  as  large  as  1  Hz  for  long  paths,  particularly  close  to  evening 
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or  daybreak.  Error  correction  coding  can  be  used  with  great  benefit  for  this  case.* 
Also,  Doppler  shift  correction  of  individual  modes  is  a  distinct  possibility,  allowing  long 
averaging  times  in  the  Rake  processor  and  a  resulting  performance  improvement  at  high 
Doppler  shifts. 

Case  (b)  applies  to  the  non-disturbed  HF  channel  also,  but  includes  the  likelihood 
that  some  mode  overlap  may  occur  under  special  propagation  conditions.  E.g.,  the 
2-hop  E  layer  may  come  close  to  the  1-hop  F2  layer  on  summer  afternoons.  For  those 
multipath  components  in  which  the  modes  overlap,  fluctuation  in  the  impulse  response 
will  occur  at  the  Doppler  difference  between  modes.  If  the  energy  in  these  fluctuating 
components  is  small  compared  to  the  non-overlapping  components,  little  performance 
loss  will  result  for  a  system  designed  as  if  Case  (a)  applied.  However,  the  use  of 
interleaving  in  addition  to  coding  would  provide  some  protection  in  the  event  that 
the  fading  energy  is  high.  Doppler  shift  correction  would  have  limited  benefit  for  the 
overlap  paths. 

Case  (c)  applies  to  the  disturbed  T1F  channel.  Interleaving  and  coding  can  provide 
benefits  here  The  statistical  model  for  the  mode  is  quite  different  from  that  of  Cases 
(a)  and  (b).  A  useful  model  here  would  be  the  complex  Gaussian  WSSUS  (wide-sense¬ 
stationary  uncorrelated  scattering)  model  [7]. 

In  the  performance  analyst  it  was  assumed  that  the  transmitter  filter  limited  the 
transmission  bandwidth  W  to  a  value  of  l/A,  where  1/A  is  the  chip  rate  and  A  is  the 
chip  duration.  This  degree  of  filtering  stretches  the  impulse  response  of  the  composite 
channel  filter  slightly,  but  does  not  degrade  the  performance  of  the  system  due  to 
the  Rake  processor.  With  this  bandwidth  limitation,  it  is  possible  to  represent  the 
composite  channel  by  means  of  a  tapped  delay  line  with  taps  spaced  A  apart  [7].  This 
representation  is  shown  in  figure  1.  The  transmitted  signal  is  represented  as  an  impulse 
train  with  areas  aqmq  transmitted  at  the  rate  1/A,  where  aq  is  the  chip  modulation 
and  mq  is  the  data  modulation.  In  this  paper  we  confine  our  attention  to  binary  phase 
shift  keying  (BPSK)  data  and  chip  modulation.  The  *-th  complex  weight  in  figure  1  is 
given  by 

bt{t)  =  h{t,iA)  (2-1) 

where  h(t ,  £)  is  the  time  variant  impulse  response  of  the  composite  channel.  A  maximum 
delay  of  M  A  is  shown  (M+l  taps)  to  accommodate  the  multipath  spread  of  the  channel. 
The  model  in  figure  1  is  applicable  to  Cases  (a),  (b),  and  (c)  discussed  above,  the 
difference  appearing  in  the  behavior  of  the  complex  tap  weights  {6,(f)}.  At  the  channel 
output  a  sampler  is  shown  representing  a  (complex)  A/D  conversion  operation  at  a  rate 
of  1/A. 


*A  phenomenon  called  twinklinq  bv  Booker  et  a]  [5]  [R]  can  occur  in  the  transition  Case?  (a)  and 
(b)  to  Case  (c).  Then  the  shallow  fading  becomes  modified  by  the  introduction  of  a  superimposed  slow  random 
modulation  with  occasional  deep  fades.  Interleaving  may  make  the  use  of  coding  effective  during  twinkling. 
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Figure  1 .  Formulation  of  Modem  Input  Using  Tapped  Delay  Line  Model  of  Composite  Channel 


SECTION  3 
RAKE  MODEMS 


The  transmitted  signal  consists  of  a  sequence  of  spread  spectrum  data  symbols 
containing  ±1  (BPSK)  modulation  according  to  the  output  of  the  binary  encoder.  Each 
symbol  is  random  appearing,  being  spread  in  bandwidth  with  ±1  chip  module: i ion  at 
the  rate  1/A.  If  we  ignore  intersymbol  interference  and  assume  that  the  fading  is 
slow  compared  to  the  sum  of  the  multipath  spread  of  the  channel  and  the  data  symbol 
time  duration,  a  matched  filter  for  the  received  symbol  followed  by  a  sampler  provides  a 
minimum  error  probability  demodulator  for  white  Gaussian  noise  (e.g.,  see  [8] ).  Ignoring 
intersymbol  interference  will  cause  little  error  since  the  spread  spectrum  processing  gain 
reduces  this  interference  to  negligible  proportions.  The  slow  fading  assumption  will  be 
assumed  valid  in  this  analysis.  It  is  not  known  whether  the  residual  noise  at  the 
excisor  output  may  be  modeled  as  Gaussian.  However,  in  any  case,  the  matched  filter 
maximizes  the  output  SNR  even  when  the  noise  is  non-Gaussian  but  white. 

The  Rake  technique  attempts  to  construct  a  filter  matched  to  each  received  symbol. 
This  matched  filter  can  be  decomposed  into  the  cascade  of  a  filter  matched  to  the 
transmitted  symbol  and  a  filter  matched  to  the  channel.  In  addition,  the  matched 
filter/sampler  operation  can  be  re- interpreted  as  a  correlation  operation.  Moreover, 
various  signal  processing  elements  may  be  rearranged  to  simplify  processing.  Because 
of  these  variations,  there  are  a  large  number  of  versions  of  the  Rake  technique.  In 
this  paper  we  examine  that  particular  subclass  of  Rake  techniques  which  appears  in 
the  form  of  a  multipath  diversity  combiner.  This  structure  can  be  shown  to  be  a 
rearranged  version  of  a  Rake  modem  consisting  of  the  cascade  of  tapped  delay  line 
channel-matching  operation  with  a  correlator-implemented  signal-matching  operation. 
The  four  techniques  examined  differ  in  the  way  the  channel  is  measured  to  provide 
weights  for  the  diversity  combiner. 

Figure  2  presents  a  block  diagram  of  the  decision-directed  coherent  Rake  modem. 
In  this  version*  the  (complex)  sampled  input  process  is  passed  through  a  (shift-register) 
tapped  delay  line  with  R  incremental  delays  of  A.  Then  the  input  plus  the  R  delayed 
versions  of  the  input  are  correlated  with  a  locally-generated,  unmodulated  PN  sequence. 
The  correlation  operation  involves  conjugate  multiplication  followed  by  an  integrate  and 
dump  (I&D)  over  successive  groups  of  N  —  Tj A  samples,  where  T  is  the  transmitted 
symbol  duration.  Each  I&D  output  is  multiplied  by  a  complex  tap  weight,  and  the  results 
summed  to  produce  the  Rake  complex  combiner  output  (i.e.,  the  channel  matched  filter 
output). 

‘Alternate  versions  exist  in  which  the  complex  digital  input  and  the  local  pseudo- noise  (PN)  sequence 
exchange  positions  or  in  which  the  input  signal  and  reference  are  both  delayed.  See  [9]  for  another  version  of  a 
decision-directed  modem  in  which  the  signal  matching  is  done  with  a  matched  filter  rather  than  a  correlator. 
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Figure  2  A  Decision-Directed  Coherent  Rake  Processor.  No  Delay  Compensation 


The  real  part  is  taken  for  BPSK  demodulation.  This  output  is  fed  to  a  hard  decision 
device  and  also  to  a  deinterleaver  and  decoder.  The  hard  decisions  are  used  in  the  tap 
weight  processing  as  discussed  below. 

The  tap  weight  measurement  circuit  shown  in  figure  2  is  a  decision-directed  pro¬ 
cedure  in  which  an  attempt  is  made  to  cancel  the  data  modulation  on  the  I&D  output 
through  feedback  of  hard  data  decisions  (with  a  delay  to  compensate  for  an  assumed 
decision  delay  of  T).  If  the  data  modulation  were  cancelled,  the  noise- free  component 
left  would  be  proportional  to  the  optimal  weight  for  predetection  diversity  combining. 
Averaging  with  a  long  time  constant  digital  filter  can  reduce  the  noise  level  sufficiently 
to  make  the  combining  effective.  The  tap  weight  filter  output  is  shown  multiplied  by  a 
factor  ep  equal  to  0  or  1.  This  factor  is  used  to  eliminate  tap  contributions  which  con¬ 
tain  such  weak  signals  that  they  will  degrade  output  SNR.  The  processing  to  determine 
ep  is  not  shown. 

Because  of  the  existence  of  a  group  delay  in  the  tap  weight  filter,  plus  the  delay  of 
T  shown  in  the  tap  weight  branch,  Doppler  shift  produces  a  spurious  phase  shift  in  the 
tap  weight  that  does  not  exist  on  the  signal  at  the  I&D  output.  This  phase  error  will 
produce  an  SNR  degradation  at  the  coherent  detector  output.  If  a  compensating  delay 
could  be  introduced  in  the  signal  path  ahead  of  the  tap  weight  multiplication,  the  same 
spurious  phase  shift  would  be  introduced  on  the  signal  component.  Then  the  weighting 
process  will  remove  the  phase  shift  effect  since  the  I&D  output  is  multiplied  by  the 
complex  conjugate  of  the  tap  weight.  Figure  3  shows  a  method  for  introducing  such  a 
delay  compensation  in  a  decision-directed  coherent  Rake  modem.  Separate  weighting 
and  combining  circuits  are  necessary  for  decision-directed  data  cancellation  and  delay- 
compensated  data  combining. 

In  the  parallel- probe  Rake  modem,  both  a  non-data  modulated  probe  DSPN  se¬ 
quence  and  a  statistically-independent,  data-modulated  DSPN  sequence  are  transmit¬ 
ted.  The  received  DSPN  probe  is  used  to  estimate  the  tap  weights  for  the  Rake  com¬ 
biner.  Figure  4  presents  a  block  diagram  of  a  parallel-probe  Rake  demodulator  proces¬ 
sor.*  Two  local  PN  generators  are  used  to  correlate  against  the  received  signal  at  each 
tap  of  the  delay  line.  The  probe  correlator  and  the  data  correlator  produce  outputs 
at  the  i&D  rate.  A  digital  filter  averages  the  probe  I&D  output  to  produce  the  tap 
weight  as  showr  .  figure  4.  Note  that  because  of  the  absence  of  the  decision-directed 
operation,  the  >t.  >pensating  delay  for  minimizing  the  Doppler  shift  degradation  can  be 
inserted  betvvt  .  e  data  I&D  and  the  tap  weighting. 

The  block  di-*"-  n  of  a  serial-probe  coherent  Rake  demodulator  is  shown  in  fig¬ 
ure  5.  For  t)  s  modem  every  P- th  transmitted  symbol  is  part  of  a  data  symbol  subse¬ 
quence  known  at  the  receiver.  The  tap  gain  measurement  is  conducted  only  with  this 
subsequence,  thus  avoiding  the  use  of  decision-directed  operation.  Both  the  normal  data 

*An  alternate  version  of  a  parallel-probe  Rake  modem  was  presented  in  [10]. 
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symbol  and  the  subsequence  symbol  epochs  are  obtained  from  the  known  starting  state 
of  the  local  PN  generator  in  the  same  way  they  are  obtained  from  the  PN  generator  at 
the  transmitter.  An  optional  compensating  delay  is  shown  inserted  in  the  signal  path 
to  minimize  the  phase  error  in  the  complex  tap  gain  produced  by  mean  Doppler  shift. 
Not  shown  is  the  fact  that  the  data  subsequence  must  be  inserted  after  the  coding  and 
interleaving  and  deleted  prior  to  de-interleaving  and  decoding. 

The  serial-probe  coherent  modem  requires  little  additional  hardware  compared  to 
the  decision- directed  coherent  modem.  In  contrast,  the  parallel- probe  modem  requires 
a  substantial  increase  in  hardware.  Moreover,  the  analyses  show  that  the  serial-probe 
modem  performance  is  essentially  identical  to  that  of  the  parallel-probe  modem  when 
parameters  are  set  properly. 

The  differential  phase  shift  keying  (DPSK)  Rake  modem  processing  is  shown  in 
figure  6.  For  this  version  of  the  Rake,  the  tap  weights  are  simply  the  previous  I&D 
outputs.  No  decision-directed  data  cancellation  is  necessary  since  the  coded  binary  data 
at  the  transmitter  modulates  a  binary  phase  change  (0°,180°)  onto  the  carrier.  The 
analysis  of  this  modem  may  be  seen  to  be  identical  to  a  decision-directed  Rake  modem 
in  which  there  are  no  decision  errors  and  there  is  no  tap  gain  filtering.  Because  of  the 
noisier  tap  weights,  the  DPSK  Rake  will  have  additional  SNR  degradation.  However, 
due  to  the  small  group  delay  in  the  tap  weight  path,  the  DPSK  Rake  technique  will  be 
less  susceptible  to  Doppler  shift  and  for  high  enough  Doppler  shift  will  be  superior  to 
the  other  modems. 


10 


Figure  3.  Use  of  Compensating  Delay  in  Decision-Directed  Coherent  Rake  Operation 
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Figure  4.  A  Parallel- Probe  Coherent  Rake  Processor 
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Figure  5.  A  Serial-Probe  Coherent  Rake  Processor 
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Figure  6. 


SECTION  4 
PERFORMANCE 


As  mentioned  in  section  1,  the  performance  analysis  is  confined  here  to  the  case 
of  non-overlapping,  purely  dispersive  propagation  modes  with  the  performance  condi¬ 
tioned  on  the  strengths  and  Doppler  shifts  of  the  modes.  The  analysis  includes  the  ef¬ 
fects  of  non-ideal  tap  weight  measurements.  Aside  from  the  degrading  effects  of  Doppler 
shifts,  the  tap  weights  are  degraded  by  additive  noise  and  by  multipath-produced  self 
noise.  The  correlation  operation  at  the  output  of  each  delay  line  tap  ideally  selects  the 
contribution  from  one  multipath  component  (i.e.,  tap  weight)  of  the  tapped  delay  line 
channel  model  in  figure  1  Due  to  the  finite  averaging  time  of  the  correlator,  contribu¬ 
tions  appear  from  all  other  multipath  components,  but  at  a  reduced  level,  producing  a 
self  noise. 

In  section  A. 4. 3  of  the  appendix,  an  analysis  is  conducted  to  evaluate  multipath 
self  noise.  It  is  found  that  if  the  Rake  processor  input  SNR,  pi  is  small  enough,  the  self 
noise  can  be  neglected  at  the  Rake  combiner  output  in  relation  to  the  additive  noise 
contributions.  Specifically,  it  is  shown  that  if 

pi  <  to/2 

self  noise  may  be  neglected,  where  the  parameter 

R0 

w~T,,W 

in  which  Rs  is  the  discrete  autocorrelation  function  of  the  composite  channel  sampled 
impulse  response  (see  (2-1)  and  figure  1) 

=  bp (*)  hp+si*)  =  Y1  ***(*’ PA)  -I-  sA)  .  (4-3) 

p  p 

The  parameter  to  may  also  be  expressed  in  the  frequency  domain  as 

\m/,t)\2d/]2 


where  //(/,<)  is  the  composite  channel  transfer  function  (Fourier  transform  of  h(t,£)). 
As  may  be  seen  from  (4-2),  to  is  upper-bounded  by  1.  For  a  single  purely  dispersive 


(4-4) 


(4-2) 


(4-1) 


propagation  mode,  w  will  be  very  nearly  unity  because  the  propagation  medium  will 
have  a  nearly  constant  magnitude  over  the  1-MHz  transmission  band.  Multiple  modes 
will  cause  w  to  become  smaller  than  unity.  Detailed  numerical  evaluations  of  w  have 
not  been  carried  out  as  of  the  writing  of  this  report,  but  it  does  appear,  based  upon 
preliminary  calculations,  that  w  will  not  be  very  much  less  than  unity  for  the  wide 
band  HF  channel.  Thus,  if  pi  «<  1,  it  appears  self  noise  may  be  neglected.  The 
calculations  presented  here  assume  that  inequality  (4-1)  is  satisfied  and  self  noise  may 
be  neglected. 

In  order  to  simplify  the  error  rate  calculations,  it  was  assumed  that  the  number  of 
taps  in  the  Rake  combiner  is  much  bigger  than  1  so  that  one  may  use  the  Central  Limit 
Theorem  to  justify  normally  distributed  statistics  for  the  combiner  output.  With  such 
statistics,  the  SNR  is  the  fundamental  parameter  in  the  determination  of  error  rate. 
Thus,  the  hard  decision  error  rate  q  (no  decoding)  is  given  by 

q  =  4>(y/p) 

where  p  is  the  combiner  output  SNR  and 


As  shown  in  the  appendix,  for  the  decision-directed  Rake  with  no  decision  errors 
and  for  the  other  three  Rake  techniques,  the  S’R  at  the  Rake  combiner  output  can  be 
expressed  in  the  generic  form 


(4-5) 


(Eb/Np)  A2 

+  l)  +  D 


(4-7) 


where  Eb/N0  is  the  energy/bit/noise  power  density  at  the  Rake  processor  input.  See 
table  1  for  a  list  of  these  parameters  for  the  four  Rake  techniques  and  the  definition  of 
some  fundamental  system  parameters  (Q,  Ar7',pn,t/n,7,  P,  T,  g,  and  G(f)). 

Decision  errors  produce  two  kinds  of  degrading  effects:  a  suppression  of  signal  level 
and  an  introduction  of  additional  noise.  An  approximate  analysis  of  these  effects  (see 
A. 4. 4)  yields  the  following  expression  for  combiner  output  SNR 


_ (Eb/Np)  A2 _ 

B(TJK  +  ')  +  +  0  -  2 ?)2)  D  +  a2(Eh/N0)  A 2 


(4-8) 


where 

a  =  4  ^  '  '7m(Pm  ~  Q  )  (4-9) 
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and  7 m  is  the  autocorrelation  function  of  the  discrete  impulse  response  {  am  ;  m  = 
0, 1, . . .  }  and  Pm  is  the  probability  that  a  pair  of  hard  decision  errors  occur  spaced  m 
symbols  apart.  It  appears  (4-8)  will  be  most  accurate  for  long  integration  times.  Note 
that  q,  the  hard  decision  error  rate,  is  a  function  of  p  so  that  (4-8)  is  a  nonlinear  equation 
in  p  (or  q ).  The  values  of  p  (and  q)  are  solved  first  with  no  delay  compensation  (g  =  0). 
Then  (4-8)  is  evaluated  using  the  appropriate  value  of  g  and  the  previously  calculated 
value  of  q  for  g  =  0,  to  obtain  the  combiner  output  SNR  with  delay  compensation.  The 
missing  element  in  the  above  formulation  is  Pm.  Calculations  presented  here  assumed 
independent  bit  errors.  However,  it  is  expected  that  the  decision  feedback  will  produce 
error  bursts.  As  a  result  the  calculations  will  be  optimistic. 

Tables  2-4  present  calculations  of  performance  for  the  four  Rake  processors  based 
upon  (4-7)  and  (4-8).  To  simplify  the  presentation,  the  propagation  modes  were  as¬ 
sumed  to  have  the  same  Doppler  shift.  A  constraint  length  7  convolutional  code  with 
three-bit  quantized  soft  decision  Viterbi  decoding  was  assumed  with  an  information 
rate  of  75  bits/second.  The  tables  show  the  required  value  of  E(,/N0  to  achieve  10“5 
decoded  error  rate  for  three  values  of  Rake  taps  Nt  =  100,50,20  and  six  values  of 
Doppler  shift  v  =  0.1, 0.2, 0.4, 0.6, 0.8  and  1.0  Hz.  Since  only  as  many  Rake  taps  are 
used  as  are  necessary  to  extract  the  significant  multipath  energy,  one  may  regard  the 
corresponding  multipath  spreads  as  not  exceeding  N?  psec  (assuming  A  =  1  psec). 

The  digital  tap  gain  filter  chosen  was  a  boxcar  shape.  The  memory  or  integration 
time  of  this  filter  was  adjusted  to  minimize  the  required  SNR  at  each  Doppler  shift, 
and  for  the  parallel-  and  serial-probe  cases,  it  was  necessary  to  jointly  optimize  the 
probe  power  and  the  filter  memory.  For  reference  purposes,  note  that  Ei/N0  ~  4.3 
dB  is  required  to  achieve  10-5  error  rate  for  the  codec  with  BPSK  over  the  additive 
white  Gaussian  noise  (AWGN)  channel.  In  the  case  of  the  decision-directed  modem, 
differential  encoding  and  decoding  were  assumed  to  counteract  the  possible  (0°,180°) 
ambiguity  produced  by  the  decision-directed  operation.  In  this  case  Eb/N0  ss  4.6  dB  is 
required  for  10-5  error  rate  over  the  AWGN  channel. 

It  is  interesting  to  note  that  the  performances  of  the  parallel  and  serial-probe 
modems  are  virtually  identical.  In  fact,  at  the  end  of  section  A. 5. 2  it  will  be  shown  that 
when  the  relative  power  in  the  two  probes  and  the  time  constant  of  the  two  tap  weight 
filters  are  made  the  same,  the  resulting  values  of  combiner  output  SNR  will  be  nearly 
the  same,  assuming  slow  enough  fading.  The  DPSK  Rake  modem  is  quite  inefficient 
at  the  low  Doppler  shifts  and  even  at  1  Hz  it  is  1.7  to  2.7  dB  worse  than  the  probe 
type  modems,  depending  on  the  number  of  Rake  taps.  The  decision-directed  Rake 
shows  somewhat  better  performance  than  the  probe  type  modems  (e.g.,  0.7  to  1.5  dB 
depending  upon  the  number  of  Rake  taps  and  the  Doppler  shift).  However  this  result 
is  provisional  until  the  impact  of  assuming  independent  decision  errors  is  determined. 
The  tables  show  that  the  potential  benefits  of  correcting  for  mode  Doppler  shift  are  2- 
3  dB,  with  the  greatest  benefit  associated  with  the  largest  number  of  Rake  taps.  SNR 
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improvement  can  also  be  obtained  by  applying  techniques  for  reducing  the  dispersion 
of  individual  modes  and  thus  reducing  the  required  number  of  Rake  taps  (i.e.,  1.3-2. 3 
dB  for  the  probe  cases  if  a  reduction  from  100  to  20  /xs  dispersion  were  possible).  The 
benefits  of  delay  compensation  are  greater  for  the  probe  systems  (1.0-1. 5  dB)  than  the 
decision-directed  system  (0.3-0. 7  dB).  This  is  because  of  the  increased  degradation  due 
to  decision  errors  when  the  system  attempts  to  operate  at  a  larger  value  of  q. 

The  optimum  values  of  integration  time  are  relatively  insensitive  to  the  number 
of  Rake  taps,  Nt ■  Although  the  integration  times  are  long  for  the  low  Doppler  shifts 
(seconds),  the  curves  of  Eb/N0  vs.  integration  time  have  a  broad  minima  at  low  Doppler 
shifts.  Thus,  a  substantial  decrease  in  integration  time  is  possible  without  significantly 
increasing  the  required  Eb/N0. 
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TABLE  1 

E  COMBINER  OUTPUT  SNR  FOR  FOUR  RAKE  MODEMS. 
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OPTIMIZED  FOR  EACH  SYSTEM. 


MINIMUM  REQUIRED  VALUES  FOR  Eb/No  AND  CORRESPONDING  OPTIMUM  TAP  WEIGHT  FILTER  MEMORY 
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VITERBI  DECODING. 

THE  UPPER  ENTRY  DENOTES  PARAMETERS  AND  PERFORMANCE  WITH  DELAY  COMPENSATION; 
THE  LOWER  ENTRY  INDICATES  THE  SAME  WITHOUT  DELAY  COMPENSATION. 


MINIMUM  REQUIRED  VALUES  FOR  Eb/No  AND  CORRESPONDING  OPTIMUM  TAP  WEIGHT  FILTER  MEMORY 
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THE  UPPER  ENTRY  DENOTES  PARAMETERS  AND  PERFORMANCE  WITH  DELAY  COMPENSATION; 
THE  LOWER  ENTRY  INDICATES  THE  SAME  WITHOUT  DELAY  COMPENSATION. 


TABLE  4 
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THE  UPPER  ENTRY  DENOTES  PARAMETERS  AND  PERFORMANCE  WITH  DELAY  COMPENSATION; 
THE  LOWER  ENTRY  INDICATES  THE  SAME  WITHOUT  DELAY  COMPENSATION. 
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APPENDIX  A 

PERFORMANCE  ANALYSIS 


A.l  INTRODUCTION 

This  appendix  presents  the  performance  analysis  of  the  Rake  modems  described 
in  the  body  of  this  report.  Section  A.2  formulates  the  input  to  the  modems,  including 
the  effect  of  the  propagation  medium.  Section  A. 3  characterizes  the  Rake  combiner 
output  for  the  decision- directed  Rake  modem,  neglecting  the  effect  of  decision  errors. 
The  approach  to  error  rate  evaluation  is  developed  in  section  A. 4.  This  latter  section 
includes  an  examination  of  the  impact  of  self  noise  (A.4.3)  and  the  effect  of  decision 
errors  (A. 4. 4).  Section  A. 5  considers  the  other  three  Rake  modems. 


A.2  SAMPLED  SIGNAL  AT  MODEM  INPUT 

Assuming  that  the  noise  samples  at  the  excisor  output  are  uncorrelated  and  sta¬ 
tionary  and  that  the  percentage  of  the  band  excised  is  sufficiently  small  to  model  the 
signal  output  as  proportional  to  the  sampled  signal  input,  the  channel  model,  from  the 
transmitter  to  excisor  output,  may  be  represented  as  shown  in  figure  A.l.  Without  loss 
of  generality,  the  transmitted  signal  is  conveniently  represented  by  an  impulse  train 
with  areas  aqmq  transmitted  at  the  rate  of  W  =  1  /A  per  second,  where  A  is  the  DSPN 
(direct  sequence  pseudo-noise)  signal  chip  duration,  aq  is  the  q- th  DSPN  chip  modula¬ 
tion,  and  mq  is  the  data  modulation  superimposed  on  the  q-th  DSPN  chip.  The  data 
modulation  is  at  a  much  slower  rate  than  the  chip  modulation, 

™q+kN  =  <4;  0  <  q  <  N  -  1  (A.l) 

where  <4  is  the  data  modulation  for  the  fc-th  data  symbol.  Each  symbol  spans  a  sequence 
of  N  DSPN  chips  {dq+kN,  0  <  q  <  N  —  1}.  The  time  duration  of  a  data  symbol  is 

T  =  IV  A  (A.2) 

Starting  with  an  impulse  train  input  signal,  the  received  signal  may  be  regarded 
as  being  formed  by  passing  the  input  through  four  filters  in  cascade:  a  chip  pulse¬ 
forming  filter,  a  transmitter  filter  which  attempts  to  confine  the  transmitted  signal  to 
a  bandwidth  of  W  Hz  (say  1  MHz),  an  HF  propagation  medium  “filter”,  and  a  receiver 
anti-alias  filter.  These  four  filtering  operations  may  be  combined,  for  analysis  purposes, 
into  one  composite  filter.  Since  we  assume  the  transmitted  signal  to  be  limited  to  a 


bandwidth  W  Hz,  this  composite  filter,  in  turn,  may  be  represented  by  a  tapped  delay 
line  with  taps  spaced  A  =  1  jW  (say  1  ^s)  apart  and  with  complex  weights  equal  to 
sampled  values  of  the  composite  channel  impulse  response  [A.l].  This  representation  is 
shown  in  figure  A.l  in  which  the  t-th  complex  weight 


bi(t)  =  h(t ,  iA) 


(A.3) 


where  h(t ,  £)  is  the  composite  channel  time  variant  impulse  response.  A  maximum  delay 
of  M  A  is  shown  (M  +  1  taps)  to  accommodate  the  multipath  spread  of  the  channel. 

The  tapped  delay  line  filter  output  plus  white  noise  are  passed  through  a  brickwall 
filter  of  bandwidth  W  and  then  sampled  at  a  rate  W /sec  to  generate  the  modem  input 
samples.  The  k-th  sample  is  given  by 


Wk  =  cjt  +  njt 


(A-4) 


where  cjt  is  the  signal  sample  and  njt  is  the  noise  sample. 

We  now  relate  the  sampled  signal  values  {c/}  to  the  input  chip  modulation  sequence 
{m*a*}  and  to  the  composite  channel  tap  multipliers  {bj(t)}.  Let  z(t)  represent  the 
output  of  the  brickwall  filter,  then 


z(t)  =  ^  mjtafc  sine  W(t  —  &A) 


(A.5) 


where 


sin  nx 

sine  x  =  - 

1 KX 


The  composite  channel  filter  output  w{t) 


is  then  given  by 


(A. 6) 


M 

w(t)  =  bi(t)z(t  —  ?A) 

i=0 


(A. 7) 
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Figure  A.l.  Formulation  of  Modem  Input  Using  Tapped  Delay  Line  Model  of  Composite  Channel 


Using  (A. 5)  in  (A. 7) 


M 

akbj(t)  sine W(t  —  *A  —  kA)  (A. 8) 

j= o  k 

Changing  summation  indices 

q  =  i  +  k,  (A. 9) 

equation  (A. 8)  becomes 

»«  =  E  (E  akbq-k(t)^  sinciy(i  —  qA)  (A. 10) 

9  k 

Thus,  after  sampling  at  t  =  iA,  we  find 

M 

ci  =  mk*kh-k(lb)  =  WA)  ae~k  (A. 11) 

k  k=0 

which  is  a  discrete  convolution. 

Note  that  in  figure  A.l  the  noise  samples  {n*}  are  formed  by  passing  white  noise 
through  a  brickwall  filter  of  bandwidth  W  and  sampling  the  result  at  IT/sec.  This 
construction  produces  uncorrelated  noise  samples. 

A. 3  RAKE  COMBINER  OUTPUT 

Figure  A. 2  presents  a  block  diagram  of  the  decision-directed  coherent  Rake  modem. 
In  this  version  the  (complex)  sampled  input  process  is  passed  through  a  (shift-register) 
tapped  delay  line  with  R  incremental  delays  of  A.*  Then  the  input  plus  the  R  delayed 
versions  of  the  input  are  correlated  with  a  locally-generated,  unmodulated  PN  sequence. 
The  correlation  operation  involves  conjugate  multiplication  followed  by  an  integrate  and 
dump  (I&D)  over  successive  groups  of  N  input  samples.  Each  I&D  output  is  multiplied 
by  a  complex  tap  weight  and  the  results  summed  to  produce  the  Rake  combiner  output. 
This  output  is  fed  to  a  deinterleaver  and  decoder  and  also  to  a  hard  decision  device. 
The  hard  decisions  are  used  in  the  tap  weight  processing  as  discussed  below. 

As  the  analysis  will  show,  ideally,  each  integrate  and  dump  output  suppresses  all 
multipath  components  except  one  of  the  M  +  1  shown  in  figure  A.l  and,  as  a  result, 

*  Alternate  versions  exist  in  which  the  complex  digital  input  and  the  local  PN  sequence  exchange  positions 
or  in  which  the  input  signal  and  reference  are  both  delayed. 
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produces  an  output  identical  to  that  for  a  hypothetical  channel  having  only  one  of  the 
composite  channel  tapped  delay  line  outputs.  Thus,  the  successive  I&D  outputs  for  a 
particular  correlator  contain  a  data  signal  proportional  to  one  of  the  composite  channel 
impulse  response  samples.  A  tap  weight  measurement  circuit  attempts  to  measure  this 
complex  impulse  response  sample.  Multiplication  of  each  I&D  output  by  the  conjugate 
of  the  corresponding  tap  weight  measurement  followed  by  summation  of  all  products 
then  approximates  predetection  or  coherent  combining  of  multipath  components. 

The  tap  weight  measurement  circuit  shown  in  figure  A.2  is  a  decision-directed 
procedure  in  which  an  attempt  is  made  to  cancel  the  data  modulation  on  the  I&D 
output  through  icedback  of  hard  data  decisions  (with  a  delay  to  compensate  for  an 
assumed  decision  delay  of  T).  If  the  data  modulation  were  cancelled,  the  noise-free 
component  left  would  be  proportional  to  the  desired  impulse  response  sample.  Averaging 
with  a  long  time  constant  digital  filter  can  reduce  the  noise  level  sufficiently  to  make 
the  predetection  combining  effective.  The  tap  weight  filter  output  is  shown  multiplied 
by  a  factor  ep  equal  to  0  or  1.  This  factor  is  used  to  eliminate  tap  contributions  which 
contain  such  weak  signals  that  they  will  degrade  output  SNR. 

The  block  diagrams  for  the  other  non-decision-directed  Rake  modems  differ  from 
figure  A. 2  primarily  in  the  method  used  to  generate  the  tap  weights.  In  the  serial 
probe  Rake  modem,  a  percentage  of  the  symbols  are  modulated  with  a  known  data 
sequence.  Correlation  is  carried  out  only  with  the  known  data  subsequence.  Thus, 
decision-directed  operation  is  unnecessary  as  is  the  delay,  T,  shown  in  figure  A.2  used 
to  compensate  for  decision  delays.  In  the  parallel-probe  Rake  modem,  a  separate  uncor¬ 
related  unmodulated  DSPN  signal  is  transmitted  in  parallel  with  the  data  modulated 
DSPN  sequence.  A  separate  correlation  operation  with  the  probe  sequence  provides  the 
complex  tap  gain  values  for  weighting  the  data  I&D  outputs.  The  DPSK  Rake  modem 
is  the  simplest  of  all  and  has  a  block  diagram  that  differs  from  figure  A.2  in  that  the 
decision-directed  operation  is  removed.  Then  the  previous  I&D  output  is  used  as  the 
complex  tap  weight.  However,  DPSK  modulation  is  required  at  the  transmitter  and 
DPSK  demodulation  is  produced  by  the  tap  weight  conjugate-multiply  operation.  In 
practice,  the  decision-directed  coherent  Rake  modem  would  employ  differential  encoding 
at  the  transmitter  and  differential  decoding  at  the  receiver,  since  the  decision-directed 
operation  can  produce  180°  phase  ambiguities. 

The  detailed  analysis  presented  is  for  the  decision-directed  coherent  Rake  modem, 
first  assuming  no  decision  errors  and  then,  in  section  A. 4. 4,  including  the  effect  of 
decision  errors.  Section  A. 5  extends  the  analysis  to  the  other  three  Rake  modems. 

In  this  section  we  will  develop  analytic  expressions  for  the  combiner  output  which 
include  the  nonidealities  associated  with  tap  weight  measurement  (excluding  decision 
errors)  and  the  basic  correlation  operation.  These  nonidealities  include  additive  noise, 
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Figure  A. 2.  Decision-Directed  Coherent  Rake  Modem,  Delayed  Signal  Version 


phase  errors  due  to  Doppler  shift,  and  “self  noise”,  but  do  not  include  decision  errors. 
Self  noise  is  due  to  the  residual  signal  components  at  an  I&D  output  due  to  undesired 
multipath  components.  These  cannot  be  suppressed  completely  due  to  the  random 
nature  of  the  PN  sequence  and  the  finite  integration  time  of  the  I&D. 

The  local  PN  sequence  {am}  fed  in  parallel  to  all  correlators  is  given  by 

=  am—R  (A. 12) 


Then  each  correlator  will  ideally  select  only  that  multipath  component  synchronized 
with  am-R,  the  i?A-delayed  component  from  the  channel.  We  examine  the  quantity 
ypk,  the  I&D  output  due  to  a  correlator  input  delayed  by  pA,  assuming  integration  over 
the  k- th  symbol  time  interval 


(k+l)N-l+R 


Vpk 


_I  v 

N  ^ 


N- 1 


wm-pam~R 


m—kN+R 


1  » 

—  jy  WQ+kN+R-paq-\-kN 


9=0 


This  I&D  output  has  a  signal  component  spk  and  a  noise  component  npk 

Vpk  —  spk  T  npk 

where,  using  (A. 4), 

1  N~1 

spk  =  ~*7  yi  cq+kN+R-p  aq+kN 


9=0 
N- 1 


—  Vs  * 

npk  ~  jy  2-j  n<l+kN+R-p  aq 


q+kN 


9= 0 


The  noise  correlations  are 
N- 1  N- 1 


_  l  ' _ ^  ^  _ 

n*pknrl  =  jy2  E  E  n9 


q+kN+R-p  ns+(N+R-r  aq+kN  as+(N 


q= 0  3=0 


Since,  by  hypothesis, 


nq+kN+R~p  ni+tN+R-r  ~  &iq+kN  -p)U+tN -r)  nP 


(A. 13) 

(A. 14) 

( A.15) 

(A- 16) 

(A. 17) 

(A. 18) 


a,+kN  “i +IN  =  hl+kNH>+<H)  l°l2 


(A. 19) 
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where 


it  follows  that 


f>ab 


_  r  1  ;  a 

\  0  ;  o 


=  b 
otherwise 


npk  nrl  = 


I/N‘  p  —  r,  k  —  £ 
0  ;  otherwise 


where  we  have  defined  the  excisor  output  complex  noise  strength 

W=I 

and  the  PN  sequence  strength, 

W*=1 


(A. 20) 


(A.21) 


(A. 22) 


(A. 23) 


Thus,  the  noise  components  at  the  I&D  outputs  are  uncorrelated  and  of  equal 
strength.  Since  we  assume  the  input  noise  to  be  zero  mean,  it  also  follows  that  the 
I&D  output  noises  are  zero  mean.  One  may  argue,  by  application  of  the  Central  Limit 
Theorem,  that  the  noises  are  very  nearly  complex  Gaussian  variables  when  N  is  large. 

Consider  now  the  use  of  the  sampled  received  signal  representation  (A. 11)  in 
(A. 15), 


J  N~  1  M 

spk  =  -jy  X/  br^q  +  kN  +  R  -  p])mq+kN+R_p_T  aq+kN+R_p_r  aq+kN  (A. 24) 

5=0  r=0 

By  separating  out  the  term  for  which  r  =  R  —  p,  we  obtain 

1  N~l 

spk  =  ^  ^  bR— j,(A[<7  T  kN  R  —  p])  friq+kN  T  xpk  (A. 25) 


5=0 


Then,  upon  using  (A.l)  in  (A. 25), 

N- 1 

N 


1 

Spk  =dk  T7  bR-P  ( +  kN  +  R  -  p))  +  xpk 


(A. 26) 


5=0 


where  dk  is  the  fc-th  data  symbol  complex  modulation,  and 
j  N- 1  M 

xPk  =  a?  X/  /C  br(&[q  +  kN  +  R-p})  mq+kN+R_p_raq+kN+R_p_r  aq+kN  (A. 27) 


5=0  r=0 

r?R-p 
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If  the  channel  changes  negligibly  over  a  time  duration  equal  to  the  sum  of  the  data 
symbol  and  impulse  response  durations,  then 

bR-p(A[q  +  kN  +  R-p ])  «  bft-p(JcT);  0<q<N-l  (A.28) 

and  (A. 26),  (A. 27)  simplify  to 

spk  =  <4  bR-p(kT)  +  xpk  (A. 29) 

xpk  =  ^2  bR-p+s{kT )  Rks  (A. 30) 

j^O 

where 

1  N~l 

Rks  ~  ~jy  'y  ^  mq+kN-s  aq+kN  —  s  aq+kN  (A. 31) 

q= 0 

Equation  (A. 29)  shows  that  for  slow  fading  the  signal  component  at  each  I&D 
output  for  the  k- th  symbol  time  contains  a  desired  signal  term  equal  to  the  product  of 
the  Jt-th  data  symbol  by  a  particular  impulse  response  sample,  plus  a  term  we  shall  call 
“self1  noise.  In  (A. 30)  this  term  is  shown  to  be  equal  to  a  weighted  sum  of  contributions 
proportional  to  each  of  the  other  multipath  components,  where  the  weights  are  random 
variables  as  given  in  (A. 31).  These  random  variables  are  finite  cross- correlations  be¬ 
tween  the  locally-generated  PN  sequence  and  the  transmitted  modulated  PN  sequence. 
If  the  averaging  time  were  infinite,  ( N  =  oo),  then  the  cross-correlation  would  be  zero 
and  the  self  noise  would  vanish.  Unfortunately,  N  is  finite  and  it  is  necessary  to  eval¬ 
uate  the  effect  of  this  self  noise  on  system  performance.  Assuming  N  »  1  and  using 
Central  Limit  Theorem  arguments,  one  may  show  that  the  self  noise  terms  have  real 
and  imaginary  parts  which  are  nearly  normally  distributed. 

We  analyze  now  the  complex  tap  weight  process.  Since  this  process  is  a  discrete 
filtered  version  of  the  corresponding  I&D  output,  it  will  also  be  representable  in  terms  of 
a  desired  signal  component,  a  self  noise  component,  and  an  additive  noise  component. 
In  the  case  of  decision-directed  operations  it  will  contain  an  additional  perturbation 
due  to  decision  errors  fed  back  to  the  data  cancellation  circuit.  This  effect  will  not  be 
considered  until  section  A. 4.4.  However,  one  may  determine  that  hard  decision  error 
rates  in  figure  A. 2  of  10-2  or  less  will  not  degrade  performance  perceptibly,  while  hard 
decision  error  rates  of  10-1  will  cause  significant  degradation.  Thus  it  is  necessary 
that  the  SNR  be  sufficient  to  keep  the  error  rate  at  the  hard  decision  device  in  figure 
A. 2  not  much  less  than  10-2  in  order  to  achieve  the  performance  predicted  in  sections 
A.4.1-A.4.3.  Section  A. 5  discusses  some  alternate  Rake  concepts  in  which  a  transmitted 
probe,  either  in  parallel  or  in  series  (time-multiplexed),  is  used.  These  techniques  avoid 
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the  degradations  caused  by  decision- directed  errors  at  the  expense  of  an  increase  in 
transmitter  power  to  achieve  the  same  error  rate  performance  as  a  decision-error- free 
decision-directed  modem. 

Using  (A. 29)  and  (A. 30)  in  (A. 14),  the  typical  I&D  output  for  slow  fading  can  be 
represented  by 

Vpk  =  dk  bR-p(kT)  +  xpk  +  ripk  (A. 32) 

Assuming  perfect  decision-directed  data  cancellation,  the  tap  weight  process  is  given 
by  (see  figure  A. 2) 

oo 

9pk  =  53  am  d k-m-1  Vp^k—m—l  (A. 33) 

m=0 

where  the  impulse  response  of  the  discrete  tap  weight  filter  is 

oo 

5(0=53  amd(t-mT)  (A. 34) 

m= 0 


Two  frequent  choices  of  impulse  responses  are  the  exponential  impulse  response, 

am  =  (l-a)Qm;  a<l  (A.35) 

corresponding  to  a  first-order  recursive  filter,  and  the  boxcar  filter  response, 

0  <  m  <  Np  —  1 


— 


0  ;  elsewhere 
Both  cases  have  been  normalized  to  unity  dc  gain, 


(A. 36) 


53  Qm  -  l 


(A. 37) 


m= 0 


Using  (A. 32)  in  (A.33) 

gPk  =  i>R-p(kT  -  T)  +  xpk  +  npl t  (A. 38) 

where  the  signal  component  of  the  tap  weight  is  given  by 

OO 

i>R-p(kT)  =  ^  am  bR_p(kT  —  mT)  (A. 39) 

m=0 
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and  the  self  noise  and  additive  noise  components  by 


xpk  =  Y  am  dk-m-i  xp,k-m-l  (A-4°) 

m=0 

oo 

"pk  =  Y  ^k-m-l  "p.fc-m-l  (A-41) 

m=0 

The  Rake  combiner  complex  output  for  the  k- th  symbol  is  given  by 

R 

g;k  yrt  (A.42) 

P= 0 

where,  as  noted  previously,  we  have  assumed  that  the  Rake  combiner  has  R  + 1  taps  and 
ev  is  1  or  0  depending  upon  whether  or  not  the  tap  output  is  included  in  the  summation. 
The  term  ep  is  included  in  the  analysis  to  model  the  operation  of  “thresholding”  to 
eliminate  taps  containing  little  signal  in  order  to  improve  the  combiner  output  SNR. 
A  number  of  algorithms  are  possible  for  determining  ep.  For  example,  ep  could  be 
determined  by  comparing  a  long-term  average  over  k  of  \gpk\2  with  a  threshold  and 
setting  ep  =  0  when  this  average  is  below  the  threshold.  Using  (A. 38),  (A. 14),  and 
(A. 29)  in  (A.42), 

R 

Ck  =  Y  ep(i>*R-p(kT  -T)  +  x*k  +  hpk)  ( dkbR-.p{kT )  +  xpk  +  npk )  (A.43) 

p= o 

For  coherent  demodulation  and  BPSK  data  and  chip  modulation,  only  the  real 
part  of  Ck  is  used.  This  real  part  can  be  separated  into  a  signal  and  a  noise  part 

Re(Ck)  =  Sk  +  Nk  (A- 44) 

where 

R 

St  =  dt  Re{  Y,  tpi>R-p(W  ~  T)bR.p(kT)}  (A.45) 

P= 0 

R 

A Tk  =  Rz{  eP  i^pkXpk  +  hpknpk  +  b*R-p(kT  -  T)(xpk  +  npk) 

P= 0 

+  xpkn*pk  +  npkx'pk  +  dkbR-p(kT)(x*pk  +  ^Jt))  }  (A-46) 


35 


and  dk  =  ±1. 

When  the  number  of  terms  in  the  summation  in  (A. 46)  is  sufficiently  large,  the  sum 
of  the  real  parts  of  the  noise  product  terms  may  be  assumed  normally  distributed  due 
to  the  application  of  the  Central  Limit  Theorem.  We  have  already  pointed  out  that  the 
noise  and  self  noise  terms  may  be  modeled  as  normally  distributed.  Thus,  for  purposes 
of  evaluating  error  rates,  we  can  model  Nk  as  a  normally  distributed  random  variable. 
Since  the  noise  and  self  noise  terms  are  zero  mean  and  uncorrelated,  Nk  will  also  be 
zero  mean.  Then  the  basic  parameter  for  determination  of  error  rate  performance  is  the 

SNR  2 

Pk  =  it  (A.47) 

Nk 

Assuming  that  the  input  noise  is  stationary,  the  variation  of  pk  with  k  will  be  due  to  time 
variations  in  the  channel.  One  must  consider,  then,  whether  error  rate  predictions  should 
be  on  a  quasi-stationary  basis  where  the  SNR  is  fixed  for  the  error  rate  computation 
and  allowed  to  vary  afterward,  or  whether  the  predictions  should  include  averages  over 
the  channel  fluctuations. 


A. 4  ERROR  RATE  EVALUATION  FOR  DECISION-DIRECTED 
MODEM 

In  this  section  we  concentrate  on  the  performance  evaluation  for  the  decision  di¬ 
rected  Rake  modem.  In  sections  A.4.1-A.4.3,  we  assume  no  decision  errors  and,  for  much 
of  the  analysis,  BPSK  data  and  chip  modulation.  Section  A. 4. 4  presents  an  approximate 
analysis  of  the  effect  of  decision  errors. 


A. 4.1  Ideal  Operation 

Consider  first  the  ideal  situation  in  which  the  tap  weight  measurement  is  perfect 
and  self  noise  is  negligible,  i.e., 


bR-p(kT  -T)  =  bR_p(kT ) 

(A. 48) 

%pk  =  Xpk  —  flpk  —  0 

(A. 49) 

R 

Sk  =  JkY,  eF  Ih-p(kT)!2 

(A. 50) 

p=0 
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Nk  =  Re{  £  eP  «pJk}  (A. 51) 

p= 0 

and,  using  (A. 21)  in  calculating  N%, 

2N 

«  =  =J-  £  e*->>  |6p(*r)|2  =  **  (A. 52) 

p=0 


where  we  have  used  the  following  identity  for  a  complex  Gaussian  variable,  z, 

1  kF  (A. 53) 

and  p0k  has  been  defined  as  the  value  of  pk  for  perfect  channel  measurement. 

The  input  SNR  in  the  bandwidth  W  is  (see  equation  (A. 41)) 


Pik  = 


?  ktl2 
\  l”*l2 


Since  (see  (A. 11)) 

lc*l2  =  E  E  bk-Pbk-q 

P  9 

we  see  from  (A. 19),  (A. 20),  and  (A. 23),  that 


M 

=  £  IWM-)!2 

p=0 

Thus,  using  (A. 23), 

1  M 

*»  =  ,£  iwwii2 

p=0 

so  that 

R 

£  |Ap(*r)i2 

P'jk  =  Pit  2A  - <  ptk  2 ,V 

£  IW«')I2 

P=0 


(A. 54) 


(A. 55) 


(A. 56) 


(A. 57) 


(A. 58) 
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where  equality  occurs  on  the  extreme  right  when  ep  =  1  and  R  =  M,  i.e.,  when  no 
Rake  tap  outputs  are  excised  and  the  number  of  Rake  taps  equals  the  number  of  taps 
on  the  composite  channel  model. 

Error  rate  performance  curves  are  frequently  presented  as  a  function  of  Ef,/N0 , 
the  ratio  of  energy/bit  to  noise  power  density,  or,  equivalently,  the  SNR,  where  the 
noise  is  measured  in  a  bandwidth  equal  to  the  bit  rate.  For  a  rate  1/2  code  and  binary 
PSK  transmissions,  the  bandwidth  equal  to  the  bit  rate  is  1/2T  where  T  is  the  symbol 
duration.  The  effective  real  (one-sided)  noise  power  density  at  the  excisor  output  is 


N0  = 


2W 


so  that  EijN0  in  the  vicinity  of  the  &-th  I&D  is  given  by 

M 


(D^TElM-,1 

P=  0 


(A. 59) 


(A. 60) 


Thus,  for  the  rate  1/2  coded  binary  PSK  system,  the  combiner  output  SNR,  pol t, 
equals  E),/N0  when  the  channel  measurement  is  perfect,  all  paths  of  significant  strength 
are  combined,  and  self  noise  is  negligible.  This  is  exactly  the  SNR  that  would  exist  at 
the  BPSK  demodulator  output  for  the  AWGN  (additive  white  Gaussian  noise)  channel. 
Consequently,  the  ideal  Rake  modem  performance  with  coding  is  identical  to  that  for  the 
AWGN  channel,  independent  of  the  channel  impulse  response  shape  for  a  given  value  of 
Ef,/N0.  Of  course  for  Ef,/N0  to  be  fixed,  the  energy  in  the  impulse  response,  ^  \br(kT)\2, 
must  remain  fixed.  For  non-overlapping,  purely  dispersive  propagation  modes,  this  en¬ 
ergy  will  fluctuate  slowly  enough  that  a  quasi-stationary  performance  analysis  should 
be  adequate.  An  experimental  determination  of  the  statistics  of  \bv(kT)\2  and  the 
noise  power  I  will  then  allow  a  calculation  of  the  statistics  of  (Eb/N0)i;  and  a  prediction 
of  the  cumulative  distribution  of  minimum  error  rate  achievable  with  the  coded  Rake 
system. 

In  practice,  the  combiner  output  SNR  will  be  less  than  ( Eb/Na)k  due  to  tap  exci¬ 
sion,  noisy  and  mismatched  tap  weights  (e.g.,  due  to  Doppler  shift),  and  possibly  self 
noise.  An  increase  in  (Eb/N0)k  is  required  to  bring  the  combiner  output  SNR  to  a  value 
that  will  achieve  the  same  error  rate  as  the  ideal  system.  Assuming  that  the  combiner 
output  SNR  can  be  determined  as  a  function  of  the  input  Eb/Na ,  the  quasi-stationary 
approach  can  be  extended  to  the  case  of  imperfect  combining.  Thus,  if  p(-)  relates  the 
combiner  output  SNR  to  the  input  value  of  Eb/N0  and  p(-)  relates  the  error  rate  of 
the  rate  1/2  coded  AWGN  channel  to  Eb/N0,  the  quasi-stationary  error  rate  of  th-. 
Rake  system  is  p(p(Eb/N0)).  In  estimating  the  cumulative  distribution  of  p,  however, 
we  need  to  consider  all  the  channel  and  Rake  parameters  affecting  the  combiner  output 
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SNR.  For  example,  for  a  fixed  number  of  Rake  taps  and  a  fixed  integration  time  for  the 
tap  weight  filter,  a  prime  determinant  of  p  will  be  the  Doppler  shifts  associated  with 
the  non-overlapping  dispersive  modes.  Thus,  in  addition  to  collection  of  the  statistics 
of  impulse  response  energy,  and  excisor  output  noise  power,  one  must  collect  statistics 
of  mode  Doppler  shift. 

Note  that  the  availability  of  the  function  p(p(Eb/N0))  allows  one  to  estimate  the 
input  value  of  EijN0  required  to  achieve  a  desired  error  rate  for  a  fixed  value  of  impulse 
response  energy.  Thus,  a  rate  1/2,  constraint  length  7,  convolutional  code  with  3-bit 
soft-decision  Viterbi  decoding  will  produce  an  error  rate  of  10-5  for  Eb/N0  «  4.3  dB, 
assuming  PSK  transmission  over  the  AWGN  channel.  The  required  Eb/N0  for,  say,  10-5 
error  rate  with  realistic  Rake  operation  can  be  obtained  by  solving  the  equation 

p(Eb/N0)  =  IQ0'43  (A. 61) 


for  Eb/N0. 

When  the  quasi-stationary  approach  cannot  be  used,  the  simple  procedure  de¬ 
scribed  above  for  error  rate  evaluation  with  coding  cannot  be  used.  A  specific  analysis 
is  needed  which  models  the  interaction  of  the  fluctuations  of  signal  at  the  combiner 
output  and  the  de-interleaving/decoding  operation. 


A. 4. 2  Output  SNR  Neglecting  Self  Noise 

In  the  absence  of  self  noise,  the  output  noise  term  (A. 46)  is  given  by 
R 

Nk  -  Re  j  ^  ep(h*k  nph  +  b*R_p(kT  -  T )  npk  +  dkbR_p(kT)  hpjt)  j  (A. 62) 

p= o 


As  discussed  previously,  the  variables  npk  and  hpk  are  modeled  as  complex  Gaussian. 
Thus,  the  summations  of  the  last  two  terms  in  (A. 62)  are  also  complex  Gaussian.  Using 
the  Central  Limit  Theorem,  the  summation  involving  the  product  term  npknpk  will 
also  be  modelable  as  a  complex  Gaussian  variable  when  the  number  of  terms  in  the 
summation  is  large,  as  assumed  here.  Since  the  mean  squared  value  of  the  real  part  of 
a  complex  Gaussian  variable  equals  one-half  of  the  mean  absolute  value  squared  for  the 
variable,  we  may  use 


ep(Kk  nPk  +  bR-p(kT  ~  T )  V  +  dk  bR-p(kT)  flU) 


(A.63) 
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To  simplify  the  calculation  further,  we  note  that  since  odd  order  moments  of 
complex  Gaussian  variables  are  zero  and  n^,  are  uncorrelated  for  all  p  and  q,  all 
cross-product  terms  obtained  by  squaring  in  (A. 63)  will  vanish  after  averaging.  Thus, 
(A. 63)  simplifies  to 


R  R 

Nk  =  lYlep  k*i2  +  \  ep  k*i2 1  bR-p(kT  -  T)\ 2 


.  M- 1  _ 

+  oJ2eP  \bR-p(kT) I2  (A-64) 


Using  (A. 21)  and  (A.41), 


_  T  00 

M2  =vEaS 


Thus, 

-  /2 

\r2  —  _ _ 

k  ~  2N2 


oo  R  j  R 

(£  “-)(£  e0  +  2iv  £  e,|6*.,(tr-r)|2 

0  p=0  p=0 

+  57?  Ea~  £  *v\<">-v(kT)? 


(A.65) 


(A. 66) 


Using  (A. 66)  and  (A. 45)  in  {A. 47),  we  obtain  the  general  expression  for  output 
SNR*,  when  self  noise  is  neglected  as 


~  Re2{  £  ep  b%_p(kT  -  T)  bR_„(kT)) 

_ P=Q  _ 

R  R  R 

fat)  (4  £  e-  +  E  'p  l6«-p(fcr)l2)  +  E  ep 


(A.67) 


p=0  p=0 


Consider  the  case  in  which  the  non-disturbed  HF  channel  consists  of  distinct  propa¬ 
gation  modes  each  with  different  Doppler  shifts  occupying  mutually  exclusive  time  delay 
segments  of  the  composite  channel  impulse  response.  Then, 


(A. 68) 


‘For  simplicity,  we  drop  the  k  subscripts  on  various  SNRs,  from  this  point  on. 
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where  un  is  the  Doppler  shift  on  the  nih  mode  and  gn(t,  £)  is  the  impulse  response  of  the 
m-th*  mode  with  the  Doppler  shift  removed.  Since  we  are  assuming  very  slow  fading,  it 
is  sufficient  to  neglect  the  time  variation  of  gn(t,0  in  the  present  analysis  and  simplify 
notation  by  using 

M*.0  =  E  9„(0  e'2’"'”'  (A.69) 


From  the  definition  of  bp(kT )  in  (A. 3), 


bp(kT)  =  Y,  <7»(pA)  e’2*"'tT 


(A. 70) 


Using  (A. 70)  in  (A. 39),  the  signal  component  of  the  tapweight  becomes 


b„(kT  -  T)  =  J2  E  Sn(pA) 

m= 0  n 


(A.71) 


br(kT  -  T)  =  E  9„(pA)  e12'1'”**-’)1'  GK) 


(A.72) 


where 


«(/)  =  J  ame~’2*'mT 


(A.73) 


is  the  transfer  function  of  the  discrete  tap  weight  filter  (A. 34). 

Using  the  “non-overlapping”  mode  hypothesis,  (A.71)  and  (A.72)  lead  to 

R  R 

E  eP  ~b*-r(lT  -  r)  bR-p(kT)  =  E  «-!>  E  MpA) I2  „) 

p=0  p=0  n 


(A.74) 


R  R 

Y  ep  \bR-p(kT)\2  =  51  e^-p  51  f^(pA)i2  =  Y  En 

p= 0  p=0  n  n 


(A. 75) 


/? 

E  |i«-p(^-7-)|2  =  £  £  M?A)|2  |G(^)|2  =  ElG("»)l2  E«  (A-™) 

p=0  p=0  n  n 
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where 


E.  =  iff~(P^)l5 

P= 0 


(A.77) 


Then,  using  (A. 74)  through  (A. 76)  in  (A. 67) 

2N 
P  = 


j-(T,  E "  Re{G’M  «,2”'”T}): 


(E  «m)(^  E  <=,  +  E  £”)  +  E 

p=0  n  n 


(A- 78) 


Note  that  with  the  non-overlapping,  purely  dispersive  channel  the  output  SNR  is  inde¬ 
pendent  of  time  (A;). 

For  a  single  propagation  model,  (A. 78)  simplifies  further  to 


2  N 


P  = 


yRe2{G'(»)  e»"T} ( £ e*_„|9(PA)|2) 

V- 0 


(E  O  E  eR--  +  E  e«-pl»(PAl2)  +  lGMf  E 

p=0 


R 

l 

V- 0 


R 

r 

p=0 


(A. 79) 


where  we  have  dropped  the  n  subsc  ipt.  p  may  also  be  represented  in  the  form 

Re2  (G*{ v)  e>2*vT)  p0 


P  = 


(EOC  +  ^  +  IGWP 


(A. 80) 


P  o 


where  it  will  be  recalled  (see  (A. 58)  without  k  subscript)  that  po  is  the  output  SNR  for 
perfect  channel  measurement,  and 


R  R 

Nt  =  T:  efl-p  =  ep  <  R  +  1 

p= o  p=0 

is  the  total  number  of  “active”  Rake  taps. 

In  the  event  that  no  tap  excision  is  used  and  M  =  R  (see  (A. 60)) 

po  =  Eb/N0 


(A. 81) 


(A. 82) 
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and 


Re 2  (G*(v)  e>2*'/T)  Eb/N0 


C“">(w+1)  +  |GM|: 

For  the  exponential  filter  (A. 35) 

oo  oo  t 

E2  /,  x2  _2m  1-0 

am  =  (l-tt)  ^  a  =jrjTJ 


GM  =  (1  -  o)  £  a-  =  ,  ri  phr 


M=0 


and  for  the  boxcar  filter  (A. 36) 


E 


°ra  = 


G(j/)  =  e-jT^f'-1)TalnirTN-Z  «  e~j*v(MF-i)T  Sinci/NFT;  uT  «  1 

Af  sm  7T£/i 


(A. 83) 

(A. 84) 
(A. 85) 

(A. 86) 
(A. 87) 


Either  filter  will  adequately  reduce  the  noise  level  on  the  tap  weight  process.  How¬ 
ever,  note  that  the  boxcar  filter  has  a  linear  phase  characteristic  corresponding  to  a 
delay  ( Np  ~  l)T/2,  whereas  the  exponential  filter  has  a  nonlinear  phase  characteristic 
which  is  approximately  linear  only  over  a  small  enough  bandwidth.  This  latter  phase 
shift  can  be  modeled  by  a  group  delay  around  zero  frequency.  Because  of  the  existence 
of  these  delays,  plus  the  delay  of  T  shown  in  the  tap  weights  branch,  the  Doppler  shift 
produces  a  spurious  phase  shift  in  the  tap  weight  that  does  not  exist  on  the  signal  at 
the  I&D  output.  This  phase  error  will  produce  an  SNR  degradation  at  the  coherent 
detector  output.  If  a  compensating  delay  could  be  introduced  in  the  signal  path  ahead 
of  the  tap  weight  multiplication,  the  same  spurious  phase  shift  would  be  introduced 
on  the  signal  component.  Then  the  weighting  process  will  remove  the  phase  shift  effect 
since  the  I&D  output  is  multiplied  by  the  complex  conjugate  of  the  tap  weight. 

Figure  A. 3  shows  a  method  for  introducing  such  a  delay  compensation  in  a  decision- 
directed  coherent  Rake  modem.  Separate  weighting  and  combining  circuits  are  used  for 
decision-directed  data  cancellation  and  delay-compensated  data  combining.  Assuming 
the  phase  shift  of  the  tap  weight  filter  is  adequately  modeled  by  the  group  delay  at 
zero  frequency,  the  net  effect  of  delay  compensation  is  to  remove  the  phase  shift  in  the 
argument  of  the  real  parts  in  (A. 78),  yielding  the  result 

7^(E  igmi*.)1 

p= - r - = -  (A.88) 

+  £»)  +  E 

n  n 
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for  multiple  modes,  and 


|G(*)|Vo 

(E^)(^  +  i)  +  IGWI2 


For  a  single  mode,  equation  (A. 89)  becomes 


P  = 


|GMP  Eb/N0 


(A. 89) 


(A. 90) 


when  no  tap  excision  is  used. 

In  the  simplest  cases,  (A. 83)  and  (A. 90),  it  is  possible  to  solve  explicitly  for  the 
Eb/N0  required  to  achieve  a  specific  /),  e.g.,  p  =  4.3  dB.  Thus  in  the  case  of  no  delay 
compensation,  equation  (A. 83),  the  Ef,/N0  required  to  achieve  a  desired  p  is  found  to 
be 

(E«™  +  |GMI2)? 

(F,b/N0)Tt,  -  2Rc,  jG,(l/)e)2,„7  >  + 


°l  + 1 GMP)  2p2  +  HE°l)(R  +  l)pRe2 

2  Re2  {G{v)e^l,T} 


(A.9I) 


In  the  case  of  the  boxcar  filter,  use  of  (A. 86)  and  (A. 87)  in  (A. 91)  produces  the 


result 


(Eb/N0)req  — 


(■p  +  sine2  vNpT'jPo 


(2sinc2  vN fT){cos2,kv(N p  +  1  )T) 


+ 


\j iyWf  +  s'nc2  p2  +  8 (R  +  l)(p/NF)(s'mc2  vNfT )  cos2  kv(Nf  +  1  )T 


(2  sine2  vPT)(cos2,kv(P  +  1)T) 


(A. 92) 
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Figure  A.3.  Use  of  Compensating  Delay  In  Decision- Directed  Coherent  Rake  Operation 


With  group  delay  compensation 


( Eb/N0)req  = 


(Ef^+jGWI  2)p 

2  |G(i/)|2  cos2  (0(u)  —  2,kvtg) 


+ 


yfe  arn  +  l<^(^)|2)2  P2  +  8  (J2  aj)  (R  +  l)p  \G{u)\ 2  cos2  (0(|/)  -  2tutq) 

2  |G(y)|2  cos2  (0(j/)  —  2ttutg)  ^ 

where  0(i/)  is  the  phase  of  G(v)  and  tg  is  the  group  delay  of  the  tap  weight  filter  at 
zero  frequency.  As  mentioned  previously,  for  the  boxcar  filter  we  are  able  to  make  the 
argument  of  the  cosine  vanish  for  all  v,  while  for  the  exponential  filter  we  can  only 
make  the  argument  negligible  for  u  sufficiently  small.  Thus,  for  the  boxcar  filter  we  can 
write, 


(Eb/N0)Teq  = 

p  +  sine2  vNpT^jp  +  yj +  sine2  vNpT Y p2  +  8(R+  l)(p/Np)  sine2  vNpT 


2  sine2  vNpT 


(A. 94) 


A. 4. 3  Output  SNR  Including  Self  Noise 
But  Neglecting  Doppler  Shift 

The  expression  for  noise  at  the  combiner  output  after  coherent  detection,  includ¬ 
ing  self  noise,  is  given  by  equation  (A.46).  Evaluation  of  the  strength  of  the  noise  is 
considerably  more  involved  when  self  noise  is  present.  While,  with  self  noise  absent,  we 
were  able  to  argue  that  the  complex  combiner  output  was  nearly  a  complex  Gaussian 
variable  for  a  large  number  of  taps  combined,  we  may  not  always  do  so  in  the  present 
case.*  Thus,  identifying  Z *  as  the  complex  combiner  output  noise,  we  cannot  use  (A. 53) 
as  was  done  in  the  absence  of  self  noise.  Instead,  we  must  use  the  identity 

jvf  =  WzT)  =  \  Izn5  +  \  flefzj}  (A. 95) 

for  computation  of  noise  power. 

In  the  case  of  complex  Gaussian  variables 

Zf  =  0  (AM) 

‘Specifically,  when  the  DSPN  chip  modulation  and  data  modulation  are  real,  e.g.,  ±1,  a  large  number 
of  Rake  taps,  Nj,  will  only  guarantee  that  the  real  and  imaginary  parts  of  the  self  noise  are  individually 
approximatable  by  Gaussian  variables. 
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and  only  the  first  term  in  (A. 95)  need  be  computed,  as  in  the  previous  section.  Consider 
now  the  summation  representing  Z *  (see  equation  (A.46)), 

R 

Zt=EeP  fa*Pk  +  -  T )  *,k  +  hR-,(kT)dki',k 

p= o 

+  xpkKk  +  *vknpk  +  hpknpk  +  'bR-pikT  -  T)  npi  +  bR-p(kT)dkh*pk)  (A. 97) 


Formulation  of  \Zk\2  or  Zk  yields  a  large  number  of  cross-product  averages.  How¬ 
ever,  most  of  these  averages  are  zero.  It  was  already  pointed  out  in  section  A. 4. 2  that 
averages  of  odd  numbers  of  additive  noise  terms  vanish.  It  is  also  true  that  averages  of 
odd  number  of  self  noise  terms  vanish.  In  addition,  the  self  noise  terms  and  additive 
noise  terms  are  zero  mean  and  there  is  statistical  independence  between  the  self  noise 
and  additive  noise  terms.  As  a  result,  the  only  possible  non-zero  average  cross-products 
in  computation  of  |Z*|2  are  cross-products  of  the  second  and  third  terms  or  the  seventh 
and  eighth  terms  in  (A. 97).  However  we  have  already  noted  in  section  A. 4. 2  that  this 
latter  cross-product  has  a  zero  average.  We  will  also  argue  that  the  former  cross-product 
averages  may  be  neglected. 

The  self  noise  terms  xp*  and  xg*  are  each  weighted  sums  of  finite-time  averages 
or  cross-correlations  between  DSPN  sequences.  (See  (A.30)  and  (A.40),  respectively.) 
However,  while  the  duration  of  the  time  average  defining  xpk  is  N  samples  (spanning 
T  =  N A  seconds),  the  duration  of  the  time  average  defining  xg*  is,  by  hypothesis,  a 
large  multiple  of  N  samples  (spanning  a  time  interval  equal  to  the  time  constant  of 
the  tap  gain  filter).  Thus,  only  a  very  small  fraction  of  the  DSPN  chip  modulations 
entering  into  the  finite-time  averages  defining  will  be  functionally  related  to  those 
entering  into  the  finite-time  averages  defining  xp*.  Consequently,  the  correlation  coeffi¬ 
cient  between  xp*  and  x?*  must  be  quite  small  and  will  be  neglected  in  the  subsequent 
analysis. 

As  a  result  of  the  above  considerations,  each  of  the  eight  terms  in  (A. 97)  will 
be  regarded  as  zero  mean  and  uncorrelated  with  all  other  terms  for  all  values  of  p  of 
interest.  Thus, 


R 

E 

q=0 


Xpk^qk  '  J'pk'L qk 


*»k<k  +  bR-p(kT-T)bR-q(kT~T ) 


XpkXqk  +  bR-p(kT)bR-q(kT )  ■  X*kXqk 


+  M2  ■  Uy,^p,+  M2  -  fa]  +  fa)„ 


(A. 98) 


where 


("0» = 2 “™><E ^ + b  E  ‘^kT  -  r>'2 

p=0  p= o 

+  2F(£a™>Eepi6«-p(w’M 


(A. 99) 


is  the  combiner  output  noise  strength  when  self  noise  is  neglected  (right  side  of  equation 
(A. 65))  and  use  has  been  made  of  (A. 20),  (A. 21),  and  (A. 64). 


The  other  component  of 


is  somewhat  simpler, 


J(^)  =  ^e{EEePe»[W  ■  XpkXqk  +  'b*R-p(kT  -T)bR-q(kT  'T) 

L  ^  p= o  q= 0 

•  (Alw» 


We  consider  in  turn  the  evaluation  of  each  of  the  averages  in  (A. 98)  and  (A. 100). 
Using  (A. 40), 

XpkXqk  =  X/X>QmQ;n  ^k-m-l  Xp,k—m  —  l  ^k-n-1  Xq,k~n~l  (A. 101) 

m  n 

With  the  definition  of  xpk  in  (A. 30)  plus  the  definition  of  d*  in  (A.l),  (A. 101)  may  be 
expressed  in  the  form 

XpkXqk  =  Y.Y.ama"  ' 
m  n 

Y,  E  hR-P+,((k  -">-!)?)  bR-,+r((k  -  "  -  l)r)  r  (A'1"2) 

S^0  r# 0 

where  Qks  is  the  short  term  autocorrelation  of  the  transmitted  data- modulated  DSPN 
sequence 

1  N-l 

Qks  ~  TJ  X/  mq+kN-s  aq+kN-s  m*q+kN  Qq+kN  (A. 103) 

q- 0 
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For  analysis  purposes,  the  statistics  of  Qks  will  be  unchanged  if  the  data  modulation  is 
removed  from  (A. 103)  since  the  transmitted  data  stream  statistics  will  be  independent 
of  the  data  modulation.  Thus,  to  simplify  notation  we  will  use 


Qka  -  vE  A(9)ag+kN-s  a*g+kN 


(A. 104) 


where 


Mg)  ~ 


1  ;  0  <  g  <  N  -  1 

0  ;  elsewhere. 


The  required  average  in  (A.  102)  can  then  be  expressed  as 

Q'ksQlr  ~  JH  ag+kN-sa*g+kNak+eN-ra 


(A. 105) 


k+lN 


(A. 106) 


g  h 


Assuming  that  the  DSPN  chip  modulation  is  a  random,  unit  magnitude,  sequence 


ag+kN-sag+kNah+lN-rah+eN 


1  ;  r  =  s  t^O,  g+kN  ~  h+£N 
p|2;  r  =  -s  ^  0,  g+kN  =  h+£N-r  (A.107) 
0  ;  elsewhere,  except  r  ^  0,  s  ^  0 


Note  that  if  the  chip  modulation  consists  of  random  selections  from  a  discrete 
uniformly  spaced  set  of  phases,  a2  will  vanish  except  for  the  case  of  two  phases  (i.e., 
±1)  when  a2  is  unity.  Using  (A.107)  in  (A. 106) 

QLQ7r  =  %  Y,  A^)A{S  +  (*  -  ON) 


&S,  —  T  I  — o"  I  2 


a^YJM9)M9  +  {k-£)N-s) 


(A. 108) 


Since 


i  Yi  Mg)Ms  +  (*  -  ON)  =  I  *  ; 


1  ;  k  =  i 
0  ;  k  /  i 


(A. 109) 


i  £  A(g)A(g  +  (*  -  £)W  -  ,)  =  Tri  ( i—lL-M.)  (A.l  10) 
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where 


It  follows  that 


Tri(x) 


1  -  jx|  ;  |x|  <  1 

0  ;  |x|  >  1 


If  (A. 112)  is  used  in  (A. 102),  it  is  found  that 
XpkXqk  =  TF  E  <*m  b*R-p+s((k  -  m  -  l)T)  bR-q+s((k  ~  m  -  1)T) 


s/  0 


2|2 


+  'lv“  QmOCn  b*R-p+s((k  -™~  l)T)  bR~q~S((k  - 

tn  n  s/0 


m 


Tri 


.  (s  —  (n  —  m)N' 


N 


(A.lll) 


(A. 112) 


m 

(A. 113) 


In  order  to  simplify  the  evaluation  of  (A.  11 3),  we  shall  ignore  the  possible  distorting 
effects  produced  by  time  variations  of  the  composite  channel  impulse  response  that  are 
significant  over  a  time  interval  of  the  order  of  the  time  constant  of  the  tap  gain  filter. 
Such  distorting  effects  must  be  kept  to  a  minimum,  anyway,  to  have  high  performance 
even  in  the  absence  of  self  noise.  Thus,  their  omission  in  the  analysis  of  self  noise  should 
only  produce  second-order  errors  in  calculation  of  modem  performance.  The  slow  time 
variation  assumption  is  modeled  by  using 


bp((k  —  m  —  1)T)  =  bp(kT)  =  bp  (A. 114) 

in  (A. 113),  resulting  in 

Kk*7k  =  ^E«rn)  (J2b*R-P+sbR-q+s  +  |a2|2J]r(-^)6^-p+s6/?_(7_Ji)  (A. 115) 

s/0  s/0 


where 


and 


r(x) 


oo 

Tri(x  —  n) 


—  OO 

n>y 


In  = 


0 


(A. 116) 


50 


is  the  discrete  autocorrelation  function  of  the  tap  gain  filter.  Note  that  r(£)  is  a 
linearly  interpolated  and  up-sampled  (by  a  factor  of  N)  version  of  the  tap  gain  filter 
autocorrelation  function  7 n,  normalized  to  70  =  Ylam-  This  normalization  produces 

r(o)  =  1. 

Consider  now  the  evaluation  of  xpkx*k.  Using  (A. 30),  we  find 

s4=EE‘m.‘^,P!:  (a. 117) 

3^0  Tqt 0 


where  Rks  is  given  by  (A. 31).  Proceeding  as  in  the  case  of  evaluation  of  QksQfr ,  we 
find  that 

^^=^+R2^»,-.Tri(iF)  (A.118) 

Thus, 

VS  =  J  E  *«-.+.  +  l“T  4  E  Tri(^)'-«-p+,  <>«-,-*  (A.119) 

3^0  3^0 


We  will  use  (A. 115)  and  (A.119)  to  obtain  a  representation  for  the  first  term  in 
the  summation  in  (A. 98),  which  is  a  self  noisexself  noise  component  of  the  (real  part) 
combiner  output  noise  power.  To  simplify  the  notation,  we  define 

.  R  R 

SN-  x  SN  ee  -  £  Y,  e”e"  <A-12°) 

p=  0  q—0 


Then,  substituting  (A. 115)  and  (A.119)  in  (A. 120)  and  changing  the  summation  over 
p  and  q  to  simplify  notation, 


2  R  R 
p= 0  q— 0 


SN-  x  SN  =  i;+.(N+»  +  l°2|2r(i) 

s/0 

X  [£  bp+,(b-+,  +  \^ |2  Tri(i)t;_,)]  (A. 121) 


We  turn  now  to  the  first  term  in  (A.  100),  which  is  the  remaining  part  of  the  self 
noisexself  noise  term.  Analogous  to  (A. 102),  we  find 


X*pkX*qk  -  /C  X/  0tm°tn 
m  n 


^  y  ^R-p-s^R-q+r  Qk-m-l,sQ*k-n-l,T 
3/O  r^O 


(A. 122) 
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where  the  required  aveiage  in  (A. 122)  can  be  expressed  as 

QkaQ*tr  =  r.r.AtoMh)  ag+kN-sag+kNah+lN-rah+tN  (A. 123) 

9  h 

Instead  of  (A.  107),  we  have 


|a2|2;  r  =  i/  0,  g  +  kN  =  h  +  £N 

ag+kN-sa*g+kNak+tN-rah+tN  =  '  1  '»  r  = -s,  g  +  kN  =  h  +  £N  -  r  (A.124) 

0  ;  elsewhere,  except  r  ^  0,s  ^  0 


Use  of  (A.124)  in  (A. 123)  results  in 


QLQTr  =  PI2  ^  Trp  -(t  <)Af) 


(A. 125) 


Thus,  (A.  122)  can  be  represented  in  the  form 


KkKic  =  ^(E«l)(|«2rE6*-p+»6«-«+>  +  Er(5)‘MA-,-.)  (A126> 


Consider  now  the  evaluation  of 

xpkxqk  ~  EE‘«  — p+a  ^ R—q+r  RksRkr 

a/0  r^O 

Proceeding  in  an  analogous  fashion  to  the  evaluation  of  Qf;SQ{ri 

RksRkr  =  |°2|2  ’  ^  Tri(^) 


which  leads  to 


xpkxqk  —  yy 


a2!2  1  s 

=  -rr  bR-P+ s  bR-g+s  +  jfYl  Tri(Tr)  bR-p+° bR- 


(A.127) 


(A. 128) 


(A. 129) 


To  simplify  notation,  we  define  the  remaining  portion  of  the  self  noisexself  noise 
term  as 

.  R  R 

SNxSN  =  -Re{Y,T,ePe’I^.  V7t}  (A130) 

p= 0  q= 0 


52 


Substituting  (A. 126)  and  (A. 129)  in  (A. 130)  and  changing  the  summations  over  p  and 
q  to  simplify  notation, 

SNxSN  = 


2  R  R 

p=0  ?=0 


e  *<={[£  *;+«  0a2i2  *;+• + ] 

3^0 

x  [E^+»(R\+>  +  Tri(^)  &,_,)]}  (A. 131) 


j^O 


We  consider  now  the  evaluation  of  the  remaining  terms  in  (A. 98)  and  (A. 100). 
These  are  included  as  special  cases  in  our  previous  evaluations,  in  particular  (A. 115), 
(A. 119),  (A. 126),  and  (A. 129).  The  second  and  third  terms  in  (A.98)  and  (A. 100)  are 
signal  x  self  noise  terms.  We  use  the  notation 
R  R 

S  X  SN  =  -  ePeq{bR-pbR-q  XpkX*qk  +  bR-pbR-q  XpkXqk\ 

p=0  9=0 

5  x  SN  =  2  £  £  e*e*  Re{bR-PbR-q  *pl?7k  +  bR-pbR-q 
p=0  9=0 

In  (A. 132)  and  (A. 133)  we  have  used 

br(kT  -T)  =  br{kT) 

as  a  consequence  of  the  slow  channel  time  variation  assumption. 

Using  (A. 115)  and  (A. 119)  in  (A. 132)  and  changing  the  summation  over  p  and  q , 
ft  ft 

S*xSN=2 (bU*  +  lfl2|2  ] 

p=0  9=0  3/O 

+  (E  ow  [  E  bl+-  (6?+»  +  PI2  rO  V»)  ]  I  (A135) 

3/0 


(A. 132) 

(A. 133) 

(A. 134) 


while  using  (A. 126)  and  (A. 129)  in  (A. 133),  we  obtain 
ft  ft 

5x5iY=277SH  eR-peR-g  Re{b*pb*  [  £  bP+s  (|a2|%+3  +  Tri(^)V*) 


p=0  9=0 


3^0 


+ (E^jy-, [E *;+.  (i“2i2&;+. + u^k-,) } }  (A136) 


The  fourth  and  fifth  terms  in  (A. 98)  are  noisexself  noise  terms.  Using  the  definition 


T  R 

n'xSN  =  2n  |l»t|2+  (A-137) 

p=0 


we  find  from  (A. 115)  and  (A.  119)  that 

2  '  R 


N.  x  SJV=lE^slf;e«.J^6p+,(6;+J  +  p|2Tri(^)6;.<) 


p=0 


L  s/0 


+  X>;+A+.  +  l«T  r<— )*„_,) 


s/0 


(A. 138) 


A. 4. 3.1  Binary  PSK  Data  and  Chip  Modulation 

The  special  case  in  which  both  the  data  and  chip  modulation  are  binary  PSK 

results  in  _ 

|a2|2  =  1  (A. 139) 

We  shall  rewrite  the  terms  SN*xSN  and  SNxSN  using  (A. 139)  and  replacing  r(-^)  and 
Tri(^)  by 

r(^)  =  1-r4)  (A-140) 

Tri(^)  =  1  -  Tric(^)  (A.141) 

To  simplify  notation,  we  define  an  alternate  tap  excision  variable* 

fP  =  CR-p  (A. 142) 


Using  (A.139)-(A.142)  in  (A. 121), 
2  R  R 


sn-xsn  =  e  Em[E  (6;+»s.  -  r<^>  6;+a-.) 


p= 0  q= 0  s/0 


X  [J2  (bP+r  Cqr  ~  T  Mjt)  bp+rb*q-r) 


(A. 143) 


r/0 


*tn  this  alternate  definition,  the  taps  are  numbered  backwards  from  the  end  of  the  tapped  delay  line. 
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where 


(A. 144) 


Cp3  —  ^P+3  4"  ^p  — 3 

(A. 144) 

Since 

OO 

E  *»  =  E(**+ k -) 

3^0  3=1 

(A. 145) 

for  any  {As},  and 

Cpj  =  Cp,  —  s 

(A. 146) 

r(^)  =  r(^) 

(A. 147) 

Tri(^> = Tri<^r) 

(A. 148) 

We  can  convert  (A. 143)  to 

2  R  R  oo 

sn-xsn=^Y1  E  WE4.^  -  rp(^)(‘^.V.  + 

p=0  g=0  5=0 

)] 

oo 

X  [5Z  CPr  C9r  ~  Wjf)  (^P+r  ^9-r  ^P-r  ^9+r)] 

r=0 

(A. 149) 

Upon  carrying  out  the  multiplications  of  the  bracketed  terms  in  (A. 149)  and  reversing 
the  order  of  summations,  we  obtain 


SN*  xSN  = 


Trfnf  1 

N 2 


oo  oo 


2  EE  |E A<i.v 


5—1  r=l  p=0 


oo  oo  R  R 

-EE[r  C(^)  +  Tric(^)]  fpbP+sCpr  f*bg-sCir} 

s=l  r= 1  p=0  9=0 


s=l  r=l 
oo  oo 


+EEr<<v)Tri^)fle{E  fpb*p+sbp+r  f<lbq-sbq-r 


s=l  r=l 


P= 0 
R 


9=0 


+E  E  /.v.»;+r}) 

p= 0  ?=0  ' 


(A. 150) 
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An  analogous  set  of  manipulations  produces  the  representation 


r1"^  2  /  i  oo  oo  It 

SNxSN  =  ^  UEE  Re2{H  frc>pr} 

3=1  r=l  p=0 


oo  oo  K  It 

-EE[r=(|r)  +  TrWF)]  ME  /p‘;+.v  E  «-»v} 


s=l  r=l 

oo  oo  it 

Y  \  '  /•  5  \  f  r  \  r>,  (  ST^ 


+EEr^Tri^N} ME  fpb'p+*bp+r  E 

3=1  r=l  p=0  g=0 


R  R  \ 

+E  UK+.K-r  E  /.»;-.»,+,})  (a. 151) 

p=0  5=0  ' 


Upon  combining  (A. 150)  and  (A. 151),  we  obtain 


2  w  °o 

SN‘  xSN  +  SNxSN  =  10  £  £  fle2{  £  /,  C;,v} 

^  3=1  r=l  p=0 

oo  oo  D  it 

-2  EE  [r«(£> + ™c(£)]  E/p&{6^}EAfie{6Uv} 

s=l  r=l  p=0  g=0 

oo  oo  /  i? 

+ 2  E E rM>  E  *  M  w4;+.}  E a M 

3=1  r  =  l  '  p=0  5=0 

/2  R 

+ E  ^ M  v->;+.}  E  a  M  v^u})]  (A152> 


Consider  now  the  S*  x  SW  term.  Using  (A.139)-(A.142)  and  (A. 144)  in  (A. 135), 


it  it 


p=0  5=0 


V”  64E  (6P+.  C  -  Trie  (0P+,K-,)' 


+  (E^)^*;[E  (WI-  -  r‘(j)bh-b,-)])  (A-153) 


With  the  symmetries  (A.145)-(A.148)  we  can  convert  (A. 153)  to  one-sided  summations 
over  s, 


S*  xSN  = 


n  R  R  /  ou  .  v  , 

2 jv  cp*  C9*  -  Tric(^)(bp+a6g-s  +  bp-»bq+»)] 

p=0  g=0  '  4=1 

+ o>2jm;[£ - rc(^)(i.;+,v. +i;-v»)])  <A-154> 


After  inverting  the  order  of  summations  in  (A.  154), 

2  \  00  . 


5—1  J?=0 

oo 

-J  £  [(E“SS,)re(*)  +  THc(*)]H«{£/,  ^P+4  £w;-,} 

g=0 


(A. 155) 


An  analogous  development  yields 

R 


3=1  p= 0 

oo  R  R 

-4  £  f(E<4.)r«W  +  Tri.w]fie{£/,‘;w  £%-.}  (A-156) 

'  3=1  p=0  9=0 


Thus, 


5*  X  SN  +  5x  SN  =  Y,  Re2{Y,frbPcr>3} 

™  3=1  P=0 

OO  RR 

£  [(  Eam)  Tc(^)  +  Tric(^r)j  i?e{  £  fp  b*  f>p+,}  fle{£/,ip,-.} 


2_ 

yv 


5  =  1 


p=0 


(A. 157) 
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We  now  turn  our  attention  to  the  last  self  noise  term,  N*xSN.  Again  using  (A. 139)- 
(A.142)  in  (A. 138), 


N’xSN  =  E  /,  [  E  V»  4.  -  E  Jf)  K- 

p= 0  3^0 

+E^+.v-Er4)^.^ 

a^O  a/0 

Using  (A.145)-(A.148)  to  convert  (A. 158)  to  one-sided  summations  over  s, 

N*  x  SN  = 

KLal) 

Icp*I  " 

5  =  1 


(A. 158) 


N2 


K  oo  oo 

E*(E  U2-E  [r=(^)+Tric(^)]*'(v.'>;-.))  (a.159) 


p=0  3=1 

and  inverting  the  order  of  summation 
N*  xSN  = 


r/V^  2  \  00  ^  °°  & 

^(EEAW’-E  [r4>+Tri4>]  E a s-.}) 

3  =  1  p=0  3  =  1  P= 0 


(A. 160) 


A. 4.3. 2  Approximations  and  Bounds 

The  results  presented  in  section  A. 4. 3. 2  apply  to  rather  general  composite  chan¬ 
nel  impulse  responses  and  are  particularly  suitable  for  numerical  evaluations.  Here  we 
wish  to  simplify  the  results  by  making  some  reasonable  assumptions  on  HF  channel 
parameters,  SNR  parameters,  and  ideal  synchronization  and  windowing  of  the  Rake 
modem. 

The  quantities  Tri(s/jV)  and  T(s/N)  are  normalized  discrete  correlation  functions 
which  have  unity  value  at  s  =  0  and  decrease  to  zero  slowly  with  s.  Tri {s/N)  is  a 
triangle  with  a  base  of  two  symbol  durations  (the  interval  (0,N)  corresponds  to  one 
symbol  duration),  while  T(s/N)  has  a  duration  of  the  order  of  twice  the  time  constant 
of  the  tap  gain  filter,  i.e.,  many  symbol  durations.  Consider  a  typical  term  in  which 
these  factors  appear,  using  Fs  to  denote  either  T(s/N)  or  Tri(s/A) 

<?„  =  E  bp+’  6l-» F-  (A161) 

3^0 
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where  the  conjugate  signs  are  optional.  If  the  factor  b^*+a  b^*}3  drops  to  zero  over  a  range 
of  s  values  small  compared  to  the  duration  of  Fa,  then  little  error  will  be  introduced  in 
crpq  if  Fs  is  replaced  by  Fo  =  1.  The  “width”  of  b^3(bg*}3)  vs .s  varies  with  p,q  from 
zero  when  p  =  q  =  0  to  a  maximum  of  2  R  for  p  =  q  =  R. 

It  appears  that  in  most  if  not  all  cases  of  interest  the  time  constant  of  the  tap  gain 
filter  will  be  much  larger  than  the  time  spread  of  the  Rake  tapped  delay  line.  It  follows 
that  V(s/N)  will  have  a  width  much  larger  than  R  and  may  be  set  equal  to  unity  in  our 
previous  equations  with  little  error.  This  is  equivalent  to  setting  Tc(s/N)  =  0.  With 
this  approximation  the  combiner  output  noise  variance  becomes 

,  2  oo  oo  it 

n I = (n!)o +^EE  ^{E^v} 


5=1  r= 1 


q  t  >  2  oo  oo  R  R 


5=1  r=l 


i  2  00 

^EftlE/-‘;v} 


„  oo  k  n 

-  Y,  ™^)Re{Y,frb?P+-}Re{T, /.%-} 


oo  R 

EE/. 

3=1  P=0 


V  I CPS  I 


t  \  '  2  oo  R 

-  E  Tri‘<v)  EWw;-}  (A162> 


where  ( N]:)0  is  given  by  (A. 99). 

By  the  same  argument,  when  the  symbol  duration  is  sufficiently  large  compared 

to  the  Rake  window  duration,  we  may  set  Tn(s/N)  =  1  or  Tr\c(s/N )  =  0  in  the  above 

expression  to  obtain  the  further  approximation 

E2  oo  oo  R  ,  *  2 

am  \  '  \  '  n  .2  f  \  '  r  *  _  1  ,  1  +  am 


^  2  w  w  a 

i?^«).+%EE  fe2{E/. 

5=1  r= 1  p—0 


V  cpscpr  f 


oo  R 


P  \cPs\ 


3=1  p= 0 


(A. 163) 
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If  the  combiner  output  SNR  is  not  very  much  larger  than  that  necessary  to  provide 
adequate  error  rate  performance  (e.g.,  10-5),  and  the  SNR  performance  degradations 
of  the  Rake  modem  are  modest,  the  input  signal-to-noise-ratio  for  the  Rake  modem,  pi, 
will  be  small  compared  to  unity,  assuming  that  the  signal  bandwidth,  e.g.,  1  MHz,  is 
much  larger  than  the  data  symbol  rate.  Then  it  is  instructive  to  examine  the  self  noise 
power  levels  relative  to  the  ambient  additive  noise  levels  at  the  Rake  system  output. 
This  will  lead  to  a  criterion  for  the  unimportance  of  self  noise  when  pi  is  sufficiently 
small. 

From  (A. 115)  we  determine  that  the  strength  of  the  p-th  tap  gain  filter  self  noise 
is 

M*  =  ^(E  I'-k-p+.I2  + 

s^O  s^O 

2  OO  OO  J  jy 

=  TTt(  E  W2  +  R2  E  r( 

fc  =  — OO  fe— — OO 


(A. 164) 

For  binary  chip  modulation  |a2|2  =  1.  The  second  sum  vanishes  for  quaternary  phase 

shift  keying  (QPSK)  modulation  (|a2|  =  0).  The  first  term  in  the  parentheses  in 

(A.  164)  is  simply  the  energy  in  the  channel’s  impulse  response.  Ignoring  the  slow  varia¬ 
tion  of  T(-)  with  k ,  the  second  term  is  essentially  a  convolution  of  the  channel’s  impulse 
response  with  its  conjugate,  evaluated  at  a  time  2 (R  —  p).  Assuming  a  highly  dispersive 
channel,  it  will  be  quite  small  compared  to  the  first  term  and,  in  any  case,  is  upper- 
bounded  by  the  first  term.  The  last  term  will  also  be  small  compared  to  the  first  term 
in  the  highly  dispersive  case,  since  |6*|2  will  have  many  terms  of  comparable  size  and 
no  one  value  will  dominate.  Thus,  with  little  error  for  binary  PSK  chip  modulation  and 
no  error  for  QPSK  chip  modulation,  we  can  set 


\xpk\  ~ 


*2 


and,  in  any  case 


\xpk  I 


< 


AT  £  I6*!2 

^Ew’ 


(A. 165) 


(A. 166) 


From  (A.  11 9)  we  see  that  the  strength  of  the  p-th  tap  output  self  noise  is  given  by 


Fp*I 


=  ifEN- 


N\^ 

s/0 


p+sl 


+ 


I  2 |2 
\az 


i  00 

4(E 


|6fc|2  + 


l^i2 


fc  =  —  OO 


3^0 

OO 

E  Tri< 

k=—oo 
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Using  the  same  arguments  as  for  the  tap  gain  filter  output  self  noise,  with  little 
error  for  BPSK  and  no  error  for  QPSK  chip  modulation,  we  can  write 


M2  =  E  i6*!2 

(A. 168) 

and,  in  any  case, 

*P*I2  -jjY,  n2 

(A. 169) 

From  (A. 21)  and  (A. 65)  we 

know  that  the  strengths 

of  the  corresponding  noise 

terms  are 

|AP*|2  = 

(A. 170) 

K* 2  =  N 

(A.171) 

Thus,  using  the  approximations  shown, 

M2 

(A. 172) 

and,  in  any  case, 

1  A  |2 

I1”*!  <2 Pi 

M2  ~ 

(A. 173) 

where  it  will  be  recalled  that  (see 

(A.57)) 

Pi  =  7  E  l6‘l2 

(A. 174) 

Similarly, 

IXP*I  ~ 

M2 

(A. 175) 

and,  in  any  case, 

I1”1!2  <2 p, 

K*l2  ' 

(A. 176) 

Thus,  when  the  input  SNR,  pi  «  1,  the  self  noise  at  each  tap  filter  and  tap 
output  will  be  much  smaller  than  the  corresponding  additive  noise.  This  might  lead 


one  to  conclude  that  the  self  noise  can  always  be  neglected  when  pi  «  1.  However, 
this.'  is  not  necessarily  true,  since  the  combining  operation  can  build  up  the  self  noise 
more  rapidly  than  the  additive  noise,  due  to  the  correlated  nature  of  the  self  noise  from 
tap-to-tap.  However,  we  show  first  that  the  noise  xself  noise  terms  may  be  ignored  by 
comparison  to  the  noisex  noise  terms  at  the  combiner  output.  Then  we  derive  criteria 
under  which  the  other  self  noise  terms  at  the  combiner  output  may  be  neglected. 

The  complex  representation  of  the  noise  xself  noise  terms  at  the  combiner  output, 
zNxSN,  are  the  fourth  and  fifth  terms  in  (A. 97),  i.e., 

R 

ZN xSN  =  ep(Xpk  h*pk  +  Kk  npk )  (A. 177) 

p=0 

Using  the  independence  of  the  tap  weight  noise  variables  and  (A.170)-(A.176),  we 
readily  find  that 

_  R  _  _  _  _ 

\zNxSn\2  =  ep(lXP*|2  ■  K*l2  +  l^pifcl2  •  |«p*|2) 

P=  0 

R  _  _ 

~2 P,^TeP  \npk\2  ■  |npjfc|2  (A. 178) 

p=0 

The  noisex  noise  term  complex  representation  zjvx//  is 

R 

ZN*N  =  ePhpknPk  (A. 179) 

P= 0 

with  corresponding  strength, 

_  R  _  _ 

ktfx/vl2  =  eP  I h*vk\2  '  K*l2  (A. 180) 

P= 0 

Thus, 

\zn*sn\2  ~  2Pi\zn*n\2  (A. 181) 

and  the  noise  xself  noise  term  may  be  neglected  when  pi  «<  1. 

We  now  develop  criteria  under  which  the  other  self  noise  terms  may  be  neglected. 
First  we  compare  the  strengths  of  the  signalxself  nuiae  teinia,  ^Xi5Jv,  v'.TIi  i,u<j  strength 
of  the  signal x  noise  terms,  zsxN >  of  the  combiner  output. 
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The  strength  of  zgxSN  is  just  twice  the  value  of  the  previously-defined  term  S*  xSN 
as  may  be  verified  by  examination  of  (A. 132).  Thus,  from  (A. 135), 


i**<™p  4EEu  k\  E  (v»  + Kl2  k+.  bu.  ^ ) ) 

p=0  q=0 

+  (e«9wE  (*;+.  wR2w.-.r(^>)l  (a.is2) 

s/0  ■* 

By  interchanging  the  summation  variables  p  and  q  in  the  second  triple  sum,  we 


=  EEw  %  E  w  *;+. 


|^5x5Ar|2 


p— 0  g=0 


I~9l2  R  R 


+  IvT  E  E  fv  /f  %  E  V*  (  Tri(^)  +  (E  r(^))  (A. 183) 


p=0  9=0 


Then,  changing  the  order  of  summations, 


IzSxSJvl 


=  1±¥iiElEw+. 


3^0  p=0 


+xE  ( Tri(^) + s r< i)) [ e  w] [ E  /p <A-i84> 


For  purposes  of  this  analysis  we  assume  that  the  number  of  Rake  taps  R  are  large 
enough  and  synchronization  has  been  adjusted  so  that  all  significant  energy  in  the 
impulse  response  is  passed  by  the  Rake  combiner.  Then, 


E  fp  bp  bp+s  w  E  bp  bp+a  —  Rs  ~  R-a 


(A. 185) 


where  Rs  is  the  autocorrelation  function  of  the  sampled  composite  channel  impulse 
response.  Thus,  using  (A. 185), 

te^P  «  1  +-^  E  l*>!2+  IT E  (™(£)  +  (E<*)r(*))l*.P  (A. 186) 
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The  last  sum  is  less  than  the  first,  since 


Tri(^)  +  (E^)r(^)<l+Eo 


Thus, 


|zs*s»l2  <  (1  +  H2)  (1  +  ^°-^)  Y,  lfi<!2 


(A. 187) 


(A. 188) 


where  the  bound  will  be  very  tight  if  the  symbol  duration  is  much  larger  than  the  width 
of  the  impulse  response  autocorrelation  function. 

The  strength  of  the  signalxnoise  terms,  zsxNi  is  readily  found  to  be 


FSxJVl 


Ro 


(A. 189) 


p=o 


where  we  have  used  the  same  approximation  (A. 185).  Using  (A. 188)  and  (A. 189)  and 
the  definition  of  px  (A.  174), 


_  _  £  i* 

|^5x5^|2  <  (f  +  |a2|2)/>»^^° 


w 


USxJVl 


(A. 190) 


Thus,  to  neglect  the  signal  x self  noise  terms,  it  is  necessary  that 

,  1  \  |7?o|2 

Pi  «  (7-7=52)  -oc—  - 

1  +  M  Y\K’\2 


(A. 191) 


The  ratio  on  the  left  can  be  converted  to  a  frequency  domain  expression.  With  the 
aid  of  Parseval’s  Theorem,  one  may  readily  show  that 


00  1  r 

Ei«»i2  =  as  / 1 W)!4# 

—  OO  ** 

iflol  2  =  ^[/ mn\2d; 


(A. 192) 


(A. 193) 


where  H(f )  is  the  composite  channel  transfer  function.  Then  the  inequality,  (A. 191), 
becomes 

Pi  <<  - =-k  w  (A. 194) 

l  +  |a2| 
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where  the  parameter 


1  fW/  2  „  -.2 

^  r  \n(f)\*df 

W  J-W/2 

Note  (from  (A. 192)  and  (A. 193))  that  w  is  upper-bounded  by  unity.  If  the  composite 
channel  transfer  function  were  an  ideal  rectangular  bandpass  filter  with  bandwidth  B, 
then 

w  —  BjW  ;  ideal  channel  of  bandwidth  B  (A. 196) 

In  such  a  case,  by  making  B  small  enough,  the  inequality  (A. 194)  will  eventually  be 
violated,  since  the  self  noise  becomes  ever  larger  as  B  is  made  smaller. 

Assuming  that  the  terminal  equipment  filter  bandwidth  W  =  1/A  Hz,  there  is 
no  way  that  the  composite  HF  channel  can  look  like  a  bandpass  filter  of  bandwidth 
B  «  W.  In  fact,  for  a  single  propagation  mode  of  arbitrary  dispersion,  w  will  be 
very  nearly  equal  to  unity,  and  \R3\2  will  drop  to  small  values  for  s  /  0.  The  described 
behavior  is  a  consequence  of  the  fact  that  for  a  single  propagation  mode  the  HF  channel 
will  have  an  essentially  constant  transfer  function  amplitude  over  the  W  =  1  MHz 
bandwidth,  so  that  \H(f)\2  will  be  determined  entirely  by  the  cascade  of  the  terminal 
equipment  filters  and  the  chip  pulse  filter. 

For  multiple  modes  w  will  become  smaller,  but  specific  calculations  have  to  be 
carried  out.  It  is  doubtful  that  w  will  become  much  smaller  than  unity  for  typical 
multimodal  HF  channel  characteristics.  As  a  consequence,  the  condition  pi  «<  1 
should  make  the  signal X  self  noise  term  small  compared  to  the  signal  x noise  term  at  the 
combiner  output. 

It  is  also  of  interest  to  compare  the  signal  x  self  noise  power  to  the  desired  signal 
power  |R0j2.  The  result  is 


(A. 195) 


"SxSjVp  <  (l  +  |^2|2)il+  ^  (Qm)2) 


£l*. 
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N 
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(1+  a22)(l+i:(am)2)  (l  +  MNl  +  EO*")2) 


N  Aw 


T  w 


(A. 197) 


where  T  =  N A  is  the  duration  of  the  integrate-and-dump.  Since  Aw  will  not  be  much 
less  than  unity,  it  is  clear  that  the  signal  x  self  noise  strength  will  be  small  compared 
to  the  signal  strength  whenever  N  >>  1  for  the  wide  band  HF  channel.  A  similar 
conclusion  was  reached  by  Price  [A. 2]  for  analog  transmission  by  the  Rake  technique. 


We  consider  now  the  strength  of  the  self  noisexself  noise  signal,  •2'SjVxSA'i  in  relation 
to  the  strength  of  the  noisexnoise,  at  the  combiner  output.  The  strength  of 

zSNxSN  is  Just  twice  the  value  of  the  term  SN’xSN  presented  in  (A. 121).  Thus, 

EEuE  (»;+.  v. + PiV v,  r<^>) 

p= 0  9=0  s^O 

*  E (w;+. + Pf v-6;-  Tri<^))  <a-198) 

.s^O 


Consider  first  the  case  of  QPSK  chip  modulation  (|a2|2  =  0),  for  which  (A. 198) 
may  be  expressed  in  the  form 


1 = EEw.  (*,-»  -  w  -  w  (a-199» 


PStVxSAM 


p=0  9=0 


where  we  have  used 

oo 

Rq-p  = bp+s  bq+3  (A. 200) 

~oo 

The  sequence  fp  represents  a  set  of  rectangular  Rake  “windows”  placed  over  in¬ 
dividual  HF  propagation  modes.  We  define  the  autocorrelation  function  of  this  (0,1) 
window  sequence  as 

Wk  =  Yl  fPfP+k  —  W °  —  Nt  (A. 201) 

v 

where  Nj  was  defined  as  the  total  number  of  taps  included  in  all  Rake  windows. 
Expanding  (A. 198)  and  using  (A. 185)  and  (A. 201), 

ks«xs»!2  =  -  2)  +  l«.l2)  (A. 202) 

If  the  inequality  (A. 201)  is  used,  we  obtain  the  bound 

\zSNxSN  I2  <  (‘V'r  y.  +  \Ro\')  (A. 203) 

The  strength  of  the  noisexnoise  term  is  from  (A. ISC),  (A. 170),  and  (A. 171 ), 


■NxN\ 
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(A.  204) 
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Using  (A. 204)  and  (A.194)  in  (A.203) 


_  2 _ 

1-zsjVxSArl2  <  ^  \zn*.n\2  (A. 205) 

assuming  Nr  »  1.  Thus,  in  order  for  the  self  noisexself  noise  term  to  be  neglected  in 
relation  to  the  noisex noise  term,  it  is  necessary  that 

w  »  p2  (A. 206) 

A  comparison  of  (A. 206)  with  inequality  (A.194)  reveals  that  satisfaction  of  (A.194) 
implies  satisfaction  of  (A. 206)  when  p,  <  1.  Moreover,  as  already  noted,  w  is  upper- 
bounded  by  unity.  Thus,  satisfaction  of  (A.194)  also  guarantees  that  the  noisexself 
noise  terms  may  be  neglected  relative  to  the  noisex  noise  terms.  For  the  case  of  QPSK, 
satisfaction  of  the  single  inequality 

w  »  pi  (A. 207) 

is  sufficient  to  allow  self  noise  effects  to  be  neglected  in  performance  analyses. 

Consider  now  the  case  of  BPSK  chip  modulation  (|a2|  =1).  Then  (A. 198)  can  be 
expressed  as  the  sum  of  four  terms 

i*™*s*ij  =  ^  [  E  E  E  *;+.  *.+.  E  w  bUr 

3^0  r^0 

+  E  E  m  E  bU.  \+.  E  Tri^>  w  b',-r 

3^0  r^0 

+ E  E  v.  E  w  bUr  E  r<^  *;+.  v. 

r^0  3^0 

+  E  E  M-  E  Tri(77>  W  K-r  E  ^  *,-•  1  (A'208) 

r^0  s/0 

Upon  interchanging  the  order  of  summations,  we  obtain 

ks«xS«i2  =  f  E  E  iA»ri2  +  EET4lfe(^ 

*•  s/0  r^0  3jt0  r^O 

+  EZr(w)/!e<A‘'A-’-') 

3^0  r^0 

+  EE  r<^)  Re<A'-'  A'-.-r}]  (A. 209) 
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where 

=  £/,&;+,  6 p+r  (A.210) 

But  it  may  be  shown  that  the  second,  third,  and  fourth  sums  in  (A. 209)  are 
individually  upper-bounded  by  the  first  sum.*  However,  this  first  sum  is  identical  to 
the  strength  of  the  self  noisexself  noise  term  for  QPSK  chip  modulation.  It  follows  from 
(A.202)-(A.205)  that  for  BPSK  chip  modulation 

IzswxsjvI2  <  ( -  P~  )  kivxjvl2  (A.211) 

w 

assuming  Nj  >>  1.  Thus,  to  neglect  the  self  noisexself  noise  term  in  relation  to  the 
noise  x  noise  term,  it  is  necessary  that 

w  »  (2 pt)2  (A. 212) 

However,  from  (A. 194)  we  see  that,  provided  2 p,  <  1,  satisfaction  of  (A. 194)  implies  sat¬ 
isfaction  of  (A. 212).  Thus,  combining  the  results  for  BPSK  and  QPSK  chip  modulation, 
provided  2 p,  <  1,  we  can  neglect  self  noise  contributions  when 

Pi  «  - 2  w  (A. 213) 

1-fa2 


A. 4. 3. 3  Output  SNR  and  Error  Rate  Evaluation 

for  Binary  PSK  Data  and  Chip  Modulation 

We  consider  here  the  evaluation  of  modem  performance  when  self  noise  cannot 
be  neglected.  The  self  noise  terms  derived  previously  ignore  any  degradations  due  to 
Doppler  shift.  This  is  consistent  with  high  bit  rate  operation,  where  the  tap  gain  filter 
time  constant  can  be  made  small  compared  with  the  reciprocal  Doppler  shift  and  yet 
keep  the  SNR  degradation  due  to  additive  noise  acceptably  small.  We  consider  then 
the  expression  for  p,  (A. 47),  the  SNR  at  the  coherent  combiner  output,  assuming  slow 
fading  and  ignoring  the  effect  of  Doppler  shift.  Then,  using  (A. 50),  (A. 99),  and  (A.  162) 
we  may  express  p  in  the  form 


Po 


i  +  E 


al(l  +  24 l)  +  2A,r^a" 


N 


P  o 


poo 


*For  example,  use  the  results  in  Chapter  X  of  [A. 3]. 


(A. 214) 


68 


where  po  is  the  value  of  p  for  perfect  channel  measurement  (see  (A. 52)  and  (A. 58))  and 
Poo  depends  on  the  self  noise  parameters  A\,B\jC\,  given  by 
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- 2  ETri^>fie{  E  w»+>}M  E  w»-.} )  (a. 217) 
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(A. 218) 


Poo  — 


(A. 219) 


The  parameter  E  is  proportional  to  the  received  signal  power  after  tap  excision. 

Computation  of  the  Ef,/N0  required  to  achieve  a  given  error  rate  proceeds  as  fol¬ 
lows.  First,  the  necessary  value  of  po  is  obtained  by  solving  (A. 214)  for  po  given  the 
value  p  needed  to  achieve  the  desired  error  rate.  This  value  of  p  has  been  shown  to 
be  equal  to  the  value  of  Ef,/N0  required  to  achieve  the  desired  error  rate  when  the 
modem  and  codec  are  operating  over  the  AWGN  channel.  As  mentioned  previously,  for 
the  rate  1/2  constant  length  7  coder  and  Viterbi  decoder  using  three-bit  soft  decisions, 
Ef,/ N0  =  4.3  dB  is  required  for  10-5  bit  error  rate  assuming  binary  PSK.  Thus,  for 
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10  5  bit  error  rate,  one  would  solve  (A. 214)  for  po  given  p  =  4.3  dB.  The  relationship 
between  E(,/N0  and  po  is  (see  (A. 52)  and  (A. 60)), 


(A. 220) 


Thus,  (A. 220)  allows  computation  of  the  required  Ef,/N0  given  the  required  value  of  po 
to  achieve  the  desired  error  rate.  If  very  little  of  the  multipath  energy  is  excluded  by 
the  tap  excision  procedure,  the  required  Eb/N0  will  be  nearly  equal  to  the  required  p0. 
In  the  case  where  there  is  too  much  multipath  to  be  handled  by  the  number  of  Rake 
taps  available,  the  required  E{,/N0  will  exceed  po . 


Note  from  (A.214)  that  for  a  given  pattern  of  excised  taps,  the  self  noise  produces 
an  irreducible  error  rate,  since  even  when  the  transmitted  power  becomes  arbitrarily 
large  and  po  — < ►  oo,  the  coherent  combiner  output  SNR  approaches  the  limiting  value 


Poo- 


lim  p  —  poo  (A. 221) 

P0—OO 


When  poo  is  less  than  the  desired  value  of  p,  e.g.,  4.3  dB,  the  desired  error  rate  per¬ 
formance  cannot  be  reached  for  any  transmitted  power.  However,  by  increasing  the 
number  of  taps  excised,  the  value  of  poo  can  be  increased  at  the  expense  of  reducing  E, 
the  effective  received  signal  power  after  tap  excision.  There  will  be  an  optimum  number 
of  taps  to  excise  which  minimizes  the  transmitter  power  required  to  achieve  a  specified 
performance. 


A. 4. 4  Effect  of  Decision-Directed  Errors 

The  previous  analyses  assume  that  the  hard  decision  error  rate  at  the  Rake  com¬ 
biner  output  is  small  enough  to  have  negligible  effect  on  the  tap  weight  generation 
process  shown  in  figure  A. 2.  However,  the  use  of  an  efficient  codec  plus  Doppler  cor¬ 
rection  and  long  tap  weight  filter  time  constants  can  reduce  Ef,/N0  requirements  to  the 
point  that  the  hard  decision  error  rate  can  be  in  the  10-1  to  10~2  range,  even  though 
the  decoded  error  rate  is  satisfactory,  e.g.,  10~5.  The  following  analysis  is  provided  to 
assess  the  SNR  degradation  caused  by  hard  decision  errors. 

To  simplify  the  presentation,  self  noise  is  neglected  and  attention  is  confined  to 
binary  PSK  chip  and  data  modulation,  although  the  analysis  is  readily  extended  to 
include  the  more  general  cases. 


70 


Assuming  decision  errors,  the  tap  weight  for  the  p-th  tap  and  the  fc-th  I&D  output 
has  the  representation 

oo 

9pk  =  ftm  dk_m_1  yp  k-m-l  (A.222) 

m=0 

instead  of  (A. 33),  where  dk  is  the  raw  data  decision  at  the  Rake  combiner  output  after 
the  k- th  I&D  outputs  have  been  weighted  and  summed. 

Using  (A. 32)  in  (A.222)  and  ignoring  the  self  noise, 

9Pk  =  ~bR-p(kT  —  T)  +  hpk  (A.223) 

where  bR-p(kT  —  T)  is  the  signal  component  and  hpk  is  the  additive  noise  component 
of  the  tap  weight  process, 


OO 

bR-p(kT)  =  otm  h-m-1  t>R-p(kT  -  mT) 


(A. 224) 


™pk  —  ^  Otm  dk-m-\  np,fc— m— 1 
m=0 


where  the  random  sequence 

(  1  ;  no  hard  decision  error  on  k- th  symbol 

Sk  =  d*k  dk  =  \ 

(  —  1 ;  hard  decision  error  on  k- th  symbol 

The  mean  signal  is 

_  00 

h-p(kT  -T)=Y,  h-m- 1  !>R-r(kT  -  mT-T ) 

m= 0 

If  q  denotes  the  hard  decision  error  probability,  then 

6k  -  l  -  2  q 


(A. 225) 


(A. 226) 


(A. 227) 


(A. 228) 


bn-p(kT)  =  (1  -  2 q)  bR.p{kT)  (A. 229) 

where  bR-p(kT  —  T),  given  by  (A. 39),  is  the  tap  weight  signal  component  without 
decision  errors.  Thus,  one  effect  of  decision  errors  is  to  suppress  the  signal  component 


at  the  tap  weight  output  by  a  factor  (1  —  2q).  Another  is  to  introduce  a  fluctuating 
component 


bR-p(kT  -  T)  -  bR-p(kT  -  T)  =  ^  am  ejt-m-i  bR-p(kT  -  mT  -  T)  (A. 230) 

m=0 

where 

ek  =  6k~l  +  2q  (A. 231) 

Including  the  variation  of  bR-p(kT  —  mT  —  T )  with  m  over  the  time  constant  of 
the  tap  weight  filter  in  (A.230)  will  have,  at  most,  a  second-order  effect  on  our  final 
results.  For  convenience  in  the  analysis,  we  neglect  this  variation  and  use 

bR-P{kT  -  T)  -  bR-p(kT  -  T)  «  bR„p(kT  -  T)  pk  ( A.232) 

where 

OO 

pk  —  'y  ]  Otm  eJt— m  —  1  (A. 233) 

m=0 

Then, 

bR-P(kT)  =  (1  -  2 q)  bR_p(kT )  T  pk  bR_p(kT)  (A.234) 

The  Rake  combiner  complex  output  for  the  k- th  symbol  is  given  by 
R 

Ck  =  ep(b*R-p{kT  -T)  +  hpk)(dk  bR-p(kT)  +  npk )  (A.235) 

p= o 

With  coherent  demodulation  and  BPSK  modulation,  only  the  real  part  of  Ck  is  impor¬ 
tant.  Thus,  we  define 

Re{Ck}  =  Sk+Nk  (A. 236) 

where  Sk,Nk  are  the  mean  signal  and  noise  components,  respectively,  of  the  coherent 
demodulated  combiner  output.  Here  we  assign  the  fluctuating  portion  of  the  signal  due 
to  previous  data  decision  errors  as  part  of  the  combiner  output  noise.  Thus, 

R 

St  =  dkYl  eP  h-p(>=T  ~  T)  bn-p(kT) 

P= 0 

R 

=  dk(  1  -  2 q)  He{  Y,ePiR-p(kT  ~  T)bR-p(kT)}  (A.237) 

P= 0 
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R  R  R 

Nk  =  #e{  ^epnphnvk  +  (l-2q+n)^ev^-v(kT-T)nvk+dk'y^epbR_v{kT)n*pk} 

p= 0  p=0  p=0 

R 

+  dkUk  Re{  Y1  ePb*R-p(kT  ~  T )  bR-p(kT)}  (A. 238) 

p=0 


In  the  previous  analysis  without  decision  errors  it  was  argued  that  N *  could  be 
modeled  as  a  normally  distributed  random  variable  via  Central  Limit  Theorem  argu¬ 
ments  based  upon  the  assumption  that  the  number  of  non-excised  Rake  taps  is  large. 
Except  for  the  last  term  in  (A. 239)  the  same  argument  applies.  The  randomness  in  the 
last  term  is  due  entirely  to  the  random  variable  e*.  As  (A.230)  shows,  this  variable  is  a 
filtered  stream  of  ±1  hard  decision  error  variables  with  the  mean  value  1  —  2 q  removed. 
In  order  to  carry  out  a  simple  analysis,  we  shall  have  to  assume  the  filter  time  constant 
is  long  enough  and  the  error  rate  is  high  enough  for  the  Central  Limit  Theorem  to 
be  invoked.  Then  we  can  assume  e*  is  also  normally  distributed.  Computation  of  error 
rates  then  can  be  carried  out  if 


(A. 239) 


is  determined. 


Proceeding  as  in  section  A. 4. 2,  we  find 
_  _  R  j2 

Nl  =  4  Re1{'£epb‘R.r(kT  -  T)bR_r(kT)}  +  ^(Eo&JJVr 

P= 0 

R  R 

+  ((1  -  2qf  +  4)  1.  £  e„  \bR_p(kT  -  r)|2  +  ~  (E  <&)(  £  e?  \bR.p(kT)\ 2) 

n  —  fl  n— O 


(A. 240) 

In  order  to  complete  the  analysis,  it  is  necessary  to  obtain  an  evaluation  of  the  strength 
of  t it-  From  (A. 233), 

OO  OO  OO 

P*  =  2  amQ "Cfc-m-lCJt-n-l  =  51  T(m)  R(^  (A.241) 

o  n  — oo 

where  '/( m)  is  the  aperiodic  autocorrelation  function  of  the  tap  weight  filter  and  Rc(m) 
is  the  autocorrelation  function  of  the  sequence  { e t).  But 

Il({n)  =  \{P{n)  -  q2)  (A.242) 
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where  P(n )  is  the  probability  that  a  pair  of  hard  decision  errors  occur  spaced  n  symbols 
apart.  Clearly 

P{  0)  =  q  (A. 243) 

P(oo)  =  q 2  (A. 244) 


If  we  ignore  the  burst  characteristics  of  {ejt},  i.e.,  assume  Cjt  is  an  uncorrelated 
sequence, 

P(n)  =  q2  ;  n  >  0  (A. 245) 


and 


r  =  4?(1  -  q)  Yj  a 


(A.246) 


It  seems  likely  that  such  an  assumption  will  lead  to  an  underestimate  for  p2  and  thus  an 
underestimate  of  the  SNR  degradation  caused  by  hard  decision  errors.  However,  note 
that  for  a  long  tap  gain  filter  impulse  response,  U  M2  will  be  small  compared 

to  unity  and  will  have  very  little  impact  on  the  output  SNR.  In  such  a  case,  one  need 
only  include  the  (1  —  2 q)  reduction  in  tap  gain  signal  on  the  output  SNR. 

Using  (A.246),  (A. 240),  and  (A. 237)  in  (A. 239),  one  may  obtain  an  expression  for 
p  that  depends  upon  q ,  the  hard  decision  error  probability.  However,  assuming  that  the 
combiner  output  is  normally  distributed,  we  have  the  relationship 


4>{Vp[q))  =  q 


(A. 247) 


where 


(A. 248) 


Thus,  to  actually  determine  q,  it  is  necessary  to  solve  the  nonlinear  equation  (A. 247). 

We  consider  this  solution  in  detail  only  for  the  case  of  a  single  propagation  mode 
with  a  Doppler  shift  u  Hz.  From  (A.74)-(A.77),  we  determine  that,  without  delay 
compensation, 


R  R 

J2e*'bK-P('kT  ~T^lt-r(kT)}  =  Re2{c(1.)e-'2"/r}( |6„_p( i-T) |- ) ' 

p— o  p— 0 

(A.  24  9) 

R  R 

Y,  eP\bR-p{kT  -  T) I2  =  !6V)|2  5]  ePN-p(U7  )|2  ( A. 250) 

p=o  p—o 

where  G{  )  is  the  transfer  function  of  the  tap  weight  filter. 
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The  combiner  output  SNR  for  perfect  channel  measurement,  po,  is  given  by 


P  o 


2N 

I 


R 

P=0 


(A. 251) 


Using  (A. 237),  (A. 240),  (A. 246),  and  (A.249)-(A.251)  in  (A.239),  we  obtain  the 
following  expression  for  the  boxcar  filter, 


_ (1-2 g)2  Re2{G( v)  e~>2™r}p0 _ 

[(1  -  2q)2  +  7]  |G(i/)P  +  A-  +  +  7l>°  Re2{G(v)  e->2"T) 


(A. 252) 


where  Np  is  the  duration  of  the  tap  filter  impulse  response  in  numbers  of  symbols.  For 
binary  PSK  data  modulation  and  a  rate  1/2  code,  it  was  shown  that 


^■tDwbi  I2 

(A. 253) 

so  that 

Eb  E  M W )l2 

N° 

(A. 254) 

Given  a  specific  tap  excision  procedure,  the  mode  Doppler  shift,  and  the  burst  parameter 
6,  it  is  possible  to  relate  po  and  q,  or  equivalently,  Eb/N0  and  q  (using  (A. 254)),  via  the 
solution  of  the  nonlinear  equations  (A. 247)  and  (A. 252). 

Numerical  evaluation  of  the  relation  between  po  and  q  is  most  simply  obtained  as 
follows.  Invert  (A. 247)  to  obtain 


M  -  [^_1  («)]2  (A. 255) 

Then  selection  of  a  value  of  q  yields  a  value  of  p.  However,  with  q  and  p  known,  the 
solution  of  po  from  (A. 252)  involves  no  more  than  the  determination  of  the  positive 
root  of  a  quadratic  equation. 

When  delay  compensation  is  used  via  the  structure  shown  in  figure  A. 3,  only  the 
SNR  at  the  combiner  output  for  the  data  path  changes  from  (A. 252).  The  data  path 
SNR  becomes 


Pdata 


_ (1-2,)2  IGWIV, _ 

[(1  -  2,)2  +  7]  ICMP  +  7^  +  +^/>o  IGMI2 


(A. 256) 
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In  order  to  evaluate  the  performance  of  the  coded  system  with  delay  compensation, 
it  is  necessary  to  develop  a  curve  relating  Pdata  f°  Po  or  Eb/N0.  The  desired  curve  is 
obtained  by  using  the  previously-computed  relationship  between  po  and  q  in  (A. 256). 
From  this  curve,  the  value  of  Eq/N0  required  to  achieve  a  desired  value  of  pdata-,  e.g., 
4.3  dB,  may  be  determined. 


Note  that  when  the  tap  weight  integration  time  is  long,  (Np  >>  1),  p2  may  be 
set  to  zero  in  (A. 252)  and  (A. 256),  and  l/Np  may  be  neglected,  yielding  the  simplified 
expressions 

(1  -2<,)2He2{GMe->2’"'7>„ 


(1-2«)J  1 GMP  + 


2  Np 

Np  po 


P  »  1 


(A. 257) 


Pdata 


(1  -  2 qf  |CM|2  po 


(i-2  ?m»)p  + 


2Np 

Np  po 


P  »  1 


(A. 258) 


A. 5  ANALYSIS  OF  THREE  OTHER  RAKE  MODEMS 

In  this  section  we  present  performance  analyses  for  three  non-decision-directed 
Rake  modems:  the  parallel-probe,  serial-probe,  and  the  DPSK  Rake  modems.  Sections 
A.5.1-A.5.3  considers  these  modems  individually. 

A. 5.1  The  Parallel-Probe  Rake  Modem 

In  the  parallel-probe  Rake  modem  both  a  non-data-modulated  probe  DSPN  se¬ 
quence  and  a  statistically  independent,  data- modulated  DSPN  sequence  are  transmit¬ 
ted.  The  received  DSPN  probe  is  used  to  estimate  the  tap  weights  for  the  Rake  combiner. 
Figure  A. 4  presents  a  block  diagram  of  the  parallel- probe  Rake  demodulator  processing. 
Two  local  PN  generators  are  used  to  correlate  against  the  received  signal  at  each  tap 
of  the  delay  line.  The  probe  correlator  and  the  data  correlator  produce  outputs  at  the 
I&D  rate.  A  digital  filter  averages  the  probe  I&D  output  to  produce  the  tap  weight  as 
shown  in  figure  A. 4.  Note  that  because  of  the  absence  of  the  decision-directed  operation, 
the  optional  compensating  delay  may  be  inserted  between  the  data  I&D  and  the  tap 
weighting. 

Using  the  composite  channel  model,  the  transmitted  signal  s(t)  consists  of  the  sum 
of  two  impulse  trains 

s (0  =  S  m*a*  ^  -  S  ^  ~  (A. 259) 
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Figure  A.4.  A  Parallel- Probe  Coherent  Rake  Processor 


where  7  is  the  ratio  of  the  power  in  the  transmitted  probe  signal  to  the  power  in  the 
transmitted  data  signal,  and  a is  the  probe  chip  modulation.  The  other  variables 
are  the  same  as  defined  previously  for  the  decision-directed  system.  Then  the  sampled 
signal  at  the  input  to  the  Rake  demodulator  may  be  represented  as 


1  /  7  (p) 

“*  =  7r ^c‘  +  i/ri^c‘  +n‘ 


(A. 260) 


where  Cfc  is  identical  to  the  sampled  received  signal  for  the  decision-directed  case  (A.ll), 
nfc  is  the  sampled  noise,  and 


4P)  -=  5Z  a‘?-k 


(r) 


(A. 261) 


is  proportional  to  the  sampled  received  probe  signal. 

Carrying  out  the  IfcD  operations  and  using  the  slow  fading  assumptions 
(c.f.(A.12)-(A.31)),  it  is  found  that  the  data  I&D  output  may  be  represented  as 


Vvk 


~  yr+7  dkbR-^kT^  +  ^/npy  Xpk  +  +  7  X^]  +  Upk  (A. 262) 


( /7  n\ 

where  the  new  term  in  (A. 262),  ,  is  an  additional  self  noise  term  due  to  the  trans¬ 

mitter  probe 

41p)  =  Y,  bK-P+.(kT)  4?'  (A. 263) 


in  which 


N- 1 


R(dp)  _  1  Y'  nip) 

Hks  ~  yy  2-j  an+kN  an+kN-s 


n= 0 


The  probe  I&D  output  has  the  analogous  form 


where 


4*”  =  E  4«-p+-(*7’)  47’ 

35^0 


in  which 


N- 1 


oM  =  i  V  a(p) *  u(/,) 

yy  /  ^  an+kN  an+k 


+k  N-s 


n= 0 


(A. 264) 


Vpt  =  *R-p(*"C)  +  +  7T+^  ^  +  mpt  (A-265) 


(A. 266) 


(A. 267) 
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and 


(Pd) 


*%d)  =  T,l«-’’+’(k7">Rk'> 


(A. 268) 


in  which 


N- 1 


D  (pd)  1  (p)* 

ft ks  —  jy  2-J  an+kN  an+kN-s 


n= 0 


The  tap  weight  process  gpk  is  given  by  an  average  over  i.e 


(A. 269) 


5pifc 


-  £  Qm  ypl- 


(A.270) 


Thus,  we  may  represent  this  process  as 


‘pp) 


9?k  -  ]/  rr^  +  / 1 + 7  ^p)  +  v/rpy  *** } + mpk 

where  the  filtered  signal  term  6/j_p(A:T)  has  been  defined  in  (A. 39)  and 


(A. 271) 


Jpp) 

xpk 

0 

(A. 272) 

oo 

r(pd) 

xpk 

-  V  a  x{pd) 

-  /  ,  am  xPlk—m 

0 

(A. 273) 

00 

mpk 

^  ^  np<k—m 

(A. 274) 

The  Rake  combiner  complex  output  for  the  fc-th  symbol  is  then 

C*  =  t  +  fih;**  +  TIT?4™  +  *'* 


P- 0 

x  |  -J=  dkbR-p(kT)  +  -7fL=  xpk  +  +  npk 


V1  +  7 


The  remaining  analysis  neglects  self  noise.  Then  the  real  part  of  the 
output  can  be  expressed  as 

Re{ck }  =  Sk  +  Nk 


(A. 275) 

combiner 
(A. 276) 
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where  the  signal  component  is  given  by 


Sk  =  dk 


^i>R-tlkT)bR_p{kT)] 


(A. 277) 


and  the  noise  component 


where 


Nk  =  Re{Zk] 


(A. 278) 


R  .  , - 

Zk  =  $>,  (:jr~dk  bR_p(kT)m'pk  +  b*R_p(kT)npk  +  npkm'pk )  (A.279) 

is  the  complex  noise  output. 

The  mean  squared  value  of  Nk  is  found  to  be  (c.f.  (A.63)-(A.66)) 


V  =  -  ^(£  ai)  JVt  +  5*  £  e, 


/ 2 


7  / 


it 


2/V2 


1  +  7  2W  £-Q 


it 


+  T-  T  £  “5.  £  ep  iw^ni2 

p= o 


1  +7  2/V 


(A. 280) 


where  Nj  is  the  total  number  of  unexcised  taps.  Using  (A. 277)  and  (A. 280),  the  output 
o.N’R  can  be  expressed  in  the  form 


Pk 


fh  /WO  +  l) 


;<£  <£)  ( 


2Nt 


/WO  +  7) 


+  l)  +\k 


(A. 281 ) 


where  the  factors  ,ik,  \k  depend  on  the  channel  impulse  response  shape  and  the  Rake 
window  function  <p  at  the  k- tli  hUI). 


R 

rt<  ~  |  ^  <P  bR_p{kI  )  bR_p(kI  )  j 


(A.  282) 


(A. 283) 


E  i 


Xk  = 


E e*  l6«-p(fcT)!i 


and  it  will  be  recalled  that 


p=0 


(A. 284) 


is  the  combiner  output  SNR  for  ideal  noise-free  tap  weights  in  the  case  of  decision- 
directed  operation. 

Assuming  ( E(,/N0)f;  is  defined  to  include  the  energy  in  the  probe  and  data 


E  b  VR-pikT)? 

p=e 


P„k  =  (Eb/N„)k  -  (A. 285) 

E  i 

/ 

which  is  the  same  relationship  as  for  the  decision-directed  modem. 

In  terms  of  the  simplified  notation  used  above,  the  combiner  output  SNR  for  the 
decision-directed  Rake  modem  operating  without  decision  errors  is 


where 


/  ^  _  $k  Pok 

[PkUecMr.  -  — TWTJ 


n 


(A. 286) 


(A. 287) 


E  |*«-,(t7-)|: 


E  ep\'h-p(kT)\‘ 

p=0 _ 

R 


(A.  288) 
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However,  with  our  basic  slow  fading  hypothesis 

Pk  «  Pk  (A. 289) 

Xk  ~  x*  (A. 290) 

Then,  comparing  (A. 286)  and  (A.281)  reveals  that  the  parallel-probe  modem  behaves 
like  a  decision-directed  modem  with  no  decision  errors  in  which  the  value  of  Ef,jN0  has 
been  reduced  by  the  factor  1/(1  +7)  and  the  tap  gain  filter  parameter  ^  is  divided 
by  7. 

With  delay  compensation  the  factors  Pk,lk  become  modified  to 


R 

*e2{E  e-‘W kT)h~p(kT-ra)} 

h  = - —R - - -  (A.291 ) 

[E*pN-p(M’-tg)|2]2 

p=0 


R 

E  l*x-r(*T)l2 

_  p= 0 

Xk  ~  r 

E  N-p(W  -  rG)|2 

P—0 


(A.  292) 


where  tq  is  the  group  delay  of  the  tap  weight  filter. 

Consider  now  the  special  case  of  the  non-disturbed  channel  with  Q  Doppler-shifted, 
very  slowly  varying,  non-overlapping  modes  (see  (A.66)-(A.77)  for  background).  In  this 
case  the  factors  Pk,Xk  become 


Pk 


[£  EqRe\G{uq)} 
7=1  _ 

«  p 

E  EA 

7=1 


(A.  293) 


Xk 


Q 


q=l 


(A.291) 


where  G(f)  is  the  transfer  function  of  the  tap  weight  filter,  uq  is  the  Doppler  shift  of 
the  q- th  modes,  and  Eq  is  its  strength  as  it  appears  after  tap  excision 


Eq  =  eR-p  |^(pA)|5 

p=0 


(A.295) 


where  gq(t)  is  the  impulse  response  of  the  n-th  propagation  mode  with  the  Doppler 
shift  removed.  The  slow  time  variation  of  gn(t)  has  been  suppressed. 

With  delay  compensation  the  factor  changes  to 


&  = 


[E  EqRe{G{vq)el2^TG}^ 

<7=1 _ 

r  Q  G 

[EM 


(A. 296) 


and  xfc  stays  the  same.  For  small  icg, 


g=t _ 

[EE«] 

7=1 


(A. 297) 


In  the  case  of  the  boxcar  filter,  (A. 297)  becomes  exact. 
The  simplest  case  is  a  single  mode,  for  which 

h  =  R<?{G(v) } 


(A. 298) 


With  delay  compensation 


which  is  exact  for  the  boxcar  filter. 


x*  = 


(3k  ~  |C(j/)|j 


(A. 299) 


(A. 300) 


If  we  assume  that  the  Rake  window  is  expanded  to  include  all  the  multipath  energy 
(/(’  =  A/,ep  =  1),  then  p0k  =  ( Ef,/N0)k .  For  the  single  mode  and  a  boxcar  tap  weight 
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filter  using  (A. 87),  (A. 298),  (A. 299)  in  (A. 281),  the  combiner  output  SNR  becomes 


P  = 


— j — (Eb/N0)  cos2  -m'(NF  —  1)T 
1+7 


/  sin  -kvNfT  \  2 
\Np  sin  ttuT ) 


+ 


/  sin  tti/N fT  \  2 
\N p  sin  tvT ) 


(A. 301) 


where  we  have  removed  the  k  subscripts  for  simplicity.  With  delay  compensation  the 
cos2  tci/(Nf  —  1  )T  factor  in  (A. 301)  would  be  eliminated. 

There  is  an  optimum  7,  i.e.,  allocation  of  power  to  the  probe  relative  to  the  data 
for  a  given  Eb/N0.  Thus,  while  a  large  probe  strength  improves  the  quality  of  the  tap 
weight,  it  reduces  the  data  signal  SNR  at  the  data  I&D  output.  This  optimum  will  vary 
with  propagation  conditions. 


A. 5. 2  The  Serial-Probe  Rake  Modem 


The  block  diagram  of  a  serial-probe  coherent  Rake  demodulatoi  processing  is 
shown  in  figure  A. 5.  For  this  modem  a  data  symbol  subsequence  known  at  the  re¬ 
ceiver  is  inserted  with  a  period  of  P  symbols.  The  tap  gain  measurement  is  conducted 
only  with  this  subsequence,  thus  avoiding  the  use  of  decision-directed  operation.  Both 
the  normal  data  symbol  and  the  subsequence  symbol  epochs  are  obtained  from  the 
known  starting  state  of  the  local  PN  generator  in  the  same  way  they  are  obtained  from 
the  PN  generator  at  the  transmitter.  An  optimal  compensating  delay  is  shown  inserted 
in  the  signal  path  to  minimize  the  phase  error  in  the  complex  tap  gain  produced  by 
mean  Doppler  shift.  Not  shown  is  the  fact  that  the  data  subsequence  must  be  inserted 
after  the  coding  and  interleaving  and  deleted  prior  to  deinterleaving  and  decoding. 


The  serial-probe  coherent  modem  requires  little  additional  hardware  compared  to 
the  decision-directed  coherent  modem,  in  contrast  to  the  parallel-probe  modem  which 
requires  a  substantial  increase  in  hardware.  In  addition,  the  known  data  subsequence 
provides  a  means  of  direct  measurement  of  raw  error  rates.  Moreover,  performance 
evaluations  for  specifications  presented  in  section  A.l  show  essentially  identical  perfor¬ 
mances  when  parameters  in  the  two  cases  are  related  properly  (e.g.,  see  (A. 334)  and 
(A.  335)). 


To  maintain  the  same  data  rate  as  the  previous  system,  it  is  necessary  that  the 
coded  symbol  rate,  i.e.,  the  IfcD  rate,  be  increased  by  a  factor  P/{P  —  1  ).  Thus,  a  new 
rer lured  symbol  duration  I1 
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Figure  A.5.  A  Serial-Probe  Coherent  Rake  Pi 


must  be  used.  The  number  of  chips  in  a  symbol  must  also  change  if  we  keep  the  chip 
rate  the  same.  Thus,  the  new  value  for  the  number  of  chips/symbol  is 

N>  =  (A. 302) 

Another  difference  between  the  systems  is  in  the  tap  weight  filters  and  tap  weight¬ 
ing.  The  serial-probe  digital  tap  weight  filter  gets  an  input  every  PT'  seconds  and  the 
same  tap  weight  is  applied  to  the  following  (P  —  1)  data  symbols  before  it  is  updated. 
The  weight  will  be  most  “out-of-date”  for  the  last  of  these  data  symbols.  Thus,  if  we 
represent  the  complex  combiner  output  for  the  k- th  I&D  output  as 

c*  =  £  s;kypt  (A. 303) 

and  we  assume  that  the  probe  symbols  occur  at  k  =  nP,  the  combiner  outputs  con¬ 
taining  data  are 

Ck  =  Y2  9p,nPVPk  ;  (A. 304) 

k  =  nP  +  l,nP  +  2,...,nP-t-P  —  1 
n  =  an  integer 

where,  assuming  slow  fading, 

Vpk  d-k  b p—p(kT  )  -|-  Xpk  T  Tipk  ,  k  ^  nP 
yPk  -  bR-p{kT')  +  xpk  +  nvk  ;  k  =  nP 

9p,nP  —  bR_p(nPT  )  -p  ip  np  -P  hp  np 


where 


i«_p(nPT')  =  £  bR.r((n  -  m)PT’)  a, 

m= 0 
oo 

% p,nP  ~  'y  %p,nP—  mP  C*m 


m= 0 

oo 


™p,nP  —  'y  ^  np,nP— 


m  P  O'm 


(A. 305) 
(A. 306) 
(A. 307) 

(A. 308) 

(A.  309) 

(A. 310) 


m=0 


/sing  (A. 305)  and  (A. 306)  in  (A. 304)  and  fA.276),  and  ignoring  self  noise 

R 

l^l !  =  ;  (A. 31 1 ) 

p= 0 

k  =  nP  - P  1 ,  n  P  -p  2,  ...,nP-p  P  —  1 
n  =  an  integer 
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IW4I2  = 


where  N' 


where 


2(^7 p(E  °^Nt  +  5F  E  e,l Ir-MPT)? 


(A. 312) 


k  =  nP  +  l,nP  +  2, . . .  ,nP  +  P  —  1 
n  =  an  integer 

is  given  in  (A. 302).  Thus  the  combiner  output  SNR  can  be  represented  as 

PkN  p'ok 


Pk  = 


(E  0(^  +  1)  + xu 

Pok 


(A. 313) 


k  =  nP  +  l,nP  +  2, . . .  ,nP  +  P  —  1 
n  =  an  integer 


R 

te2{Y,e’-iR-r(nPT">bK-plkT">} 


PkN  = 


P= 0 


(A. 314) 


E  ep 

v- o 


k  =  nP  +  l,nP  +  2,  ...,nP  +  P—  1 
n  =  an  integer 
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(A. 315) 


R 

£e,lWnPT')!2 

p=0 

xkN  =  — -  ; 

E  « , 

P=° 

A:  =  nP  +  l,nP  +  2, . . .  ,nP  +  P  —  1 

n  —  an  integer 

where  (using  (A. 302)) 

ON  P  —  1 

gt  =  —  E  (A. 316) 

is  the  combiner  output  SNR  for  ideal  noise-free  tap  weights  for  the  serial-probe  modem 
and  {p0k)v  's  the  combiner  output  SNR  for  the  decision  directed  system  under  the 
same  ideal  conditions  in  the  vicinity  of  the  time  kT' . 

If  we  include  the  probe  and  data  power  in  computing  Eb/N0  at  a  given  time 
instant,  then  this  value  is  unchanged  for  all  the  systems  analyzed  (assuming  the  same 
transmitter  power).  Thus, 

I2  (A. 317) 

p=o 

is  the  value  of  Eb/N0  for  the  serial-probe  system  in  the  vicinity  of  the  time  t  =  kT' . 
Thus  from  (A. 285)  and  (A.316),  we  see  that  we  can  set 

R 

E  l6«-#T')l2 

P= 0 _ 

M 

E  i6«-p(H’')i2 

p= o 

The  value  of  delay  to  use  in  delay  compensation  requires  some  additional  con¬ 
sideration  in  the  case  of  the  serial-probe  modem,  because  the  group  delay  of  the  tap 
weight  process  is  periodically  time- variable  due  to  the  weight  freeze.  Thus,  it  is  probably 
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best  to  use  a  compensating  delay  equal  to  the  average  delay  of  the  tap  weight  process. 
This  is  given  by  the  group  delay  of  the  digital  filter  per  se,  plus  half  the  time  that  the 
weight  is  frozen  with  the  result  rounded  off  to  the  nearest  symbol.  For  computation 
of  error  rates,  upper  and  lower  bounds  can  be  selected  by  choosing  the  particular  set 
of  data  symbols  for  which  the  delay  compensation  is  worst  and  best,  respectively.  The 
same  bounding  approach  can  be  used  when  delay  compensation  is  not  used.  This  would 
entail  evaluation  of  /?£«  and  Xfcn  at  the  two  sets  of  k  values  {k  =  nP  +  P  —  1}  and 
{&  =  nP  +1}. 

To  summarize,  then,  without  delay  compensation,  the  upper  and  lower  bounds  on 
error  rate  may  be  obtained  by  evaluating  the  expressions  for  /?£„  and  Xkn  in  (A. 314) 
and  (A. 315),  respectively,  at  the  two  values  {k  =  nP  -f  P  —  1}  and  {k  =  nP  +  1}, 
respectively.  With  delay  compensation,  /3*n  and  Xkn  would  be  chosen  as 

R 

fo2{  E  ^n-p("PT')b^p{kT'  -  GT')} 

0tn  - - ^ -  (A. 319) 

[X>l<’K-p(W’'-Gr'>l2f 

P- 0 

R 

XkN  =  -  (A. 320) 

E  eP  \bR-p(kT'  -  GT')\2 

p= 0 

where  G  is  the  integer  closest  to  the  sum  of  the  group  delay  of  the  tap  gain  filter  in 
units  of  symbol  duration,  tq/T\  plus  (P  —  l)/2.  Upper  and  lower  bounds  on  error  rate 
can  be  obtained  by  appropriate  selection  of  k  values.  This  will  most  likely  be  either 
k  =  nP  +  1  or  nP  +  P  —  1  for  the  upper  bound,  and  k  somewhere  halfway  between 
these  values  for  the  lower  bound. 

Consider  now  the  non-disturbed  wide  band  HF  channel  with  Q  Doppler  shifted 
but  very  slowly  varying,  non-overlapping  nodes  (see  (A.66)-(A.77)  for  background).  In 
this  case 


R 


Q 


E  ^pl>R-p(nPT')bR_pikTl)  =  Yjn-(vq)p’2^kT'-'pr"<  E, 


(A. 321) 


p=0 


7=1 


k  —  nP  +  1 , . . . ,  nP  -f  P  —  1 

n  =  an  integer 
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(A. 322) 


R  Q 

EepiwfcT)i2  =  E^ 

p=0  7= 1 

R  Q 

E  =»  fbR-r(nPT')\2  =  E  l6(^)l2  E,  ( A-323) 

p=0  7=1 

where  En  is  given  by  (A. 295)  and  G(v)  is  the  transfer  function  of  the  tap  gain  filter. 
Thus,  (A. 314)  and  (A. 315)  simplify  to 

r  ^  2 

[E  EqRe2  [G*(vq)  ei‘2*l'<l(kT'-nPT')} 
fan  =  — - q -  ;  (A. 324) 

[E  E<}2 

7=1 


k  —  nP  4- 1, . . .  ,nP  +  P  —  1 
n  =  an  integer 

Q 

E  i6wi2  E* 

7=1 

Xity  = - q -  , 

E^ 

7=1 


(A. 325) 


k  =  nP  +  1, . .  - ,  riP  +  P  —  1 
n  =  an  integer 

With  delay  compensation  (A. 320),  stays  the  same  as  (A. 325),  but  faw, 

(A. 319),  becomes  modified  to 
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(A. 326) 


Pkn  = 


£  EqRe{G*{vh)  e^q(kr-Gr-nPT')  | 
q-i _ 

X>12 

7=1 


&  =  nP  +  1, . . . ,  nP  +  P  —  l 
n  =  an  integer 


The  previous  comments  about  upper-  and  lower-bounding  error  rates  apply  to 
(A. 324)  and  (A.325). 

For  a  single  mode  without  delay  compensation 


fan  =  Re2{G*{u)  eVMkT'-nPT') }  . 

k  =  nP  +  l,nP  -f  2, . . .  ,nP  +  P  —  1 
n  =  an  integer 

XkN  =  \GM\2 

and  with  delay  compensation 

f3kn  =  Re2{G*{u)t32T^kT,-GT-nPr)}  ; 


(A. 327) 


(A. 328) 


(A. 329) 


k  =  nP  +  1  ,nP  +  2, . . .  ,nP  +  P  —  1 
n  =  an  integer 

and  Xkn  is  the  same  as  (A. 328). 

The  expression  for  G(f)  in  the  case  of  a  boxcar  filter  is  given  by  (c.f.  (A. 87)) 


G{u) 


3*v{Nf-\)PT‘  sinT tvNfPT' 
N'f  sin  7T 1/ PT1 


(A. 330) 


where  Np  is  the  memory  of  the  boxcar  filter.  This  must  be  distinguished  from  the 
impulse  response  duration  of  the  tap  weight  filter  NpPT' .  Without  delay  compensation, 
the  lowest  combiner  output  SNR  (which  upper-bounds  the  error  ra*e)  is  given  by 


(A. 327),  with  k  =  nP  +  P  —  1.  If,  in  addition,  we  assume  that  the  Rake  window  includes 
all  the  multipath  energy  ( R  =  M,ep  =  1),  (see  (A. 319)),  then 


and  the  combiner  output  SNR  becomes 


P  = 


l 


2  N- 


M 


sin  tv N LPT'  \  2 


N'f  sin  it 


vPT' ) 


(A. 331) 


(A. 332) 


where  we  have  removed  subscripts  to  simplify  notation.  With  delay  compensation,  the 
cos2  in/T'(PN'F+P—2)  term  becomes  changed  to  cos2  in/T'  ((Np  —  l)P+2(k  — nP)  —  2G). 
If  we  select  G  to  minimize  the  phase  error  and  then  select  k-nP  =  1, 2, ...  P-1  to  make 
the  resulting  error  as  large  as  possible,  we  find  that  the  term  becomes  approximately 
cos2  ttuT'P  where  we  have  taken  some  liberties  in  rounding  off  which  are  unimportant 
so  long  as  vT'  «  1,  as  we  assume. 

In  comparing  the  combiner  output  SNR  expression  for  the  single  mode  serial-probe 
case,  (A. 332),  to  the  single  mode  parallel-probe  case,  (A. 301),  we  find  that  by  selecting 


NF  =  N'F(P  -l)  (A. 333) 

7  =  (A. 334) 

the  two  expressions  become  essentially  identical,  assuming  uT  <<  1.  The  same  conclu¬ 
sion  results  for  the  case  of  delay  compensation 


A.5.3  The  DPSK  Rake  Modem 

The  DPSK  Rake  modem  receiver  processing  is  shown  in  figure  A. 6.  For  this  version 
of  the  Rake,  the  tap  weights  are  simply  the  previous  I&D  outputs.  No  decision-directed 
data  cancellation  is  necessary  because  the  coded  binary  data  at  the  transmitter  is 
impressed  as  a  binary  phase  change  (0°,  180°)  on  the  carrier.  The  analysis  of  this  modem 
may  be  seen  to  be  identical  to  a  decision-directed  Rake  modem  in  which  there  are  no 
decision  errors  and  there  is  no  tap  gain  filtering. 
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Thus,  neglecting  self  noise,  the  SNR  at  the  combiner  output  is  given  by  (A. 286), 
where  now  instead  of  (A. 287)  and  (A.288)  the  terms  /3k,Xk  are  given  by 


tc 

Re 2  {£* 

P= 0 

R 


(A. 335) 


Xk  =  1 


(A. 336) 


S  =  1 


(A.  337) 


For  slow  fading  and  Doppler  shifts  small  compared  to  1/T,  /?*  will  be  essentially 
unity.  Then 

n  =  — —  (A. 338) 
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Figure  A. 6.  A  DPSK  Rake  Processor 
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ERRATA 


Change  Equation  (4-2)  to  read  as  follows: 


1*0 1 


Eliminate  the  reference  to  Equation  (A.41)  above  Equation  (A. 54),  as  follows: 
The  input  SNR  in  the  bandwidth  W  is 


Pik  = 


j  M2 

j  Kl2 


Change  Equation  (A.286)  to  read  as  follows: 


(Pk)dec.dir.  — 


PkPok 


(E 


'  Pok 


+  1 


)  +Xk 


A  replacement  for  page  8  (Figure  2)  is  attached. 


(4-2) 


(A.54) 


(A.286) 


1 


Figure  2  A  Decision-Directed.  Coherent  Rake  Processor.  No  Delay  Compensation. 


